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Abstract 

In recent three-loop calculations of massive Feynman integrals within Quantum Chromody- 
namics (QCD) and, e.g., in recent combinatorial problems the so-called generalized harmonic 
sums (in short S-sums) arise. They are characterized by rational (or real) numerator weights 
also different from ±1. In this article we explore the algorithmic and analytic properties 
of these sums systematically. We work out the Mellin and inverse Mellin transform which 
connects the sums under consideration with the associated Poincare iterated integrals, also 
called generalized harmonic polylogarithms. In this regard, we obtain explicit analytic contin- 
uations by means of asymptotic expansions of the S'-sums which started to occur frequently 
in current QCD calculations. In addition, we derive algebraic and structural relations, like 
differentiation w.r.t. the external summation index and different multi-argument relations, for 
the compactification of S-sum expressions. Finally, we calculate algebraic relations for infinite 
S'-sums, or equivalently for generalized harmonic polylogarithms evaluated at special values. 
The corresponding algorithms and relations are encoded in the computer algebra package 
HarmonicSums. 



1 Introduction 



In loop calculations of Feynman Diagrams in Quantum Field Theories special classes of multiply 
nested sums and iterated integrals appear both in intermediate steps and in the final results. 
E.g., this applies both to the massless case and diagrams characterized by a single mass scale, 
or in cases two sufficiently different mass scales appear and one may consider precise power 
expansions in their ratio. We will focus in the following on single differential quantities. They 
can be characterized either by a dimensional variable x with i £ R, < i < 1 with exponential 
growth, or an integer variable N. Both representations are equivalent and are related by a Mellin 
transform!]] 

M[f(x)](n) := f x n f(x)dx. (1.1) 



J o 

The n-dependence may occur generically for Feynman diagrams with local operator insertions 
which appear, e.g., within the light-cone expansion |133|. 1141], I3"6 l l63l IT8] in deep-inelastic scatter- 
ing. As has been shown in Refs. [27 \ \29 \ [131] these Feynman integrals can always be represented 
by multiply nested sums. However, in general the specific type of these sums is not easily pre- 
dicted. It is assumed that the classes of sums do steadily extend going to higher loops and/or 
allowing for more masses and external off-shellness for a growing number of legs. 

Systematic analyses for a large class of single-scale processes at the 2-loop level in the mass- 
less case revealed )22[ [2"1] that the results for all processes could be expressed in terms of six 
basic functions in Mellin-space only. Here one allows as well to use argument-duplication and 
differential relations. The corresponding functions are multiple harmonic sumsj^] 

For a non-negative integer n and nonzero integers a» (1 < i < k) the harmonic sums [19, 123J 
are defined as 

c sign(ai) 11 sign(a fc )* fc 

&ai,... ) o fc W - 2^ • {L - Z) 

n>h>i 2 >--->i k >l H % k 

If ai 7^ 1, the limit n — > oo exists and the values can be expressed in terms of multiple zeta 
values, resp. Euler-Zagier values [56 1 1139] ; f° r further properties and an extensive list of literature 
we refer to [28]. Likewise, more general classes of sums also result into new classes of special 
numbers. 

In physics applications harmonic sums emerge for the first time in case of the fermion-fermion 
anomalous dimension (the Mellin transform of the corresponding splitting function) [60] as the 
classical harmonic sum Si(n) (also called harmonic numbers). A first non-trivial nested harmonic 
sum appears at the 2-loop level with 5_2,i(n), [67] . 

Harmonic sums satisfy a quasi-shuffle algebra |71j and can be formulated in terms of indefinite 



nested integrals. The Mellin transform (1.1) enables one to map harmonic sums to their integral 
representation in terms of weighted harmonic polylogarithms, and the inverse transform allows 
one to transform these integral representations back to harmonic sums [105] . These properties 
give rise to the analytic continuation [2j2 El ESI EH [251 ESI [27] of harmonic sums. In particular, 
one can calculate the asymptotic expansion of harmonic sums [261 [27] ; for the special case dj > 1 
see also [UJ. Moreover, various relations between harmonic sums are derived which allows to 
compactify harmonic sum expressions coming from QCD calculations. Here in particular the 
algebraic relations [901 EU 1123] induced by the quasi-shuffle algebra and the structural relations 
|27l [26] based on the differentiation of harmonic sums contribute substantially. 

Harmonic sums play a prominent role, e.g., in the field of combinatorics or number theory. 
Here the summation package Sigma |113] based on difference field theory [771 H101 11121 HH1 1115] 



1 Here / : R — > R can be distribution valued function, cf. [138] . 

2 In the physics literature their use dates back to [67j . see also [621 187 
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contributed to discover and prove identities such as [991 11091 EQl 11144 HOlj ; as a special case these 
algorithms contain hypergeometric summation |102j . Meanwhile these algorithms are combined 
with the HarmonicSums package [6j (containing the analytic and algebraic properties summarized 
above) and lead to new packages [H E3 EH E] dealing with definite multi-sums. 

In recent calculations based on these packages new classes of functions have been arising 
in intermediate results but also in the final result of given Feynman integrals. During such 
calculation but also for further considerations of the given result it is crucial that the knowledge 
of these new function classes is explored to a similar advanced level as the harmonic sums. For 
instance, the new class of cyclotomic harmonic sums has been analyzed in detail in [5]. 

In particular, the generalized harmonic sums or so-called S-sums |90j occur in higher order 
calculations P251 El El 03] . 

These sums also may emerge in connection with cyclotomic index alphabets [5] and nested 
binomial and inverse binomial weights, cf. |11| I12j. In the context of statistics or number theory 
typical calculations can be found, e.g., in [9"9" [ 11001 198], In case of pure harmonic weights they are 
defined by 



H: 



S ai ,...,a k (xi, ■ ■ ■ ,Xk;n) — ^2 7^r'''7^*r' 0-3) 



n>ii>«2>'-->ifc>l 



with positive integers en and (usually) real numbers Xi ^ (1 < i < k). Note that restricting 



±1 gives back the class of harmonic sums (1.2), or in case that Xi G {u £ C\u n = l,n E N} 



the complex representations of the cyclotomic sums [5] . There are mutual applications of these 
functions in loop calculations, for a survey see e.g. [130]. Special cases of the S'-sums are, cf. [90] : 
the harmonic sums [1231 119] and the multiple zeta values resp. Euler-Zagier values, cf. e.g. [28] . 
the multiple polylogarithms [661 135] . two-dimensional harmonic polylogarithms including the 
letters {0, 1, —1, l—z, —z] [SUES]) the harmonic polylogarithms over the alphabet {0, 1, —1} [105J, 
the Nielsen integrals over the alphabets {0, 1}, resp. {0, —1} [96[ [79] HU [80], and the classical 
polylogarithms |84| [83] [4"5l 1140] . Furthermore [HO], the generalized hypergeometric functions 
p+iFp [751 El [1201 E3 EH [58] the Appell-functions F h2 and the Kampe De Feriet functions 
[HI [151 HH I120| E3 EH EE] can be represented in terms of S'-sums. Code-implementations of four 
general algorithms being described in [90] are given in [1281 192] 1127] . In field theoretic calculations 
the e-expansion in the dimensional parameter is of importance and needs to be solved for the 
S-sums and the analytic continuation in their parameters. For special classes the e-expansion 
has been studied in Refs. [901 US E3 ESI [731 UM El EH E] . 

In [90] it is heavily exploited that these sums satisfy a quasi-shuffie algebra. S-sums form a 
Hopf algebra [721 ESI [Ell E2J ESI EE] , as has been shown in Ref. [90], see also [52] . The quasi-shuffie 
property leads to new techniques to expand certain definite nested sums over hypergeometric 
terms in terms of S-sums. In addition, it has been worked out in [35] that the S-sums at infinity 
can be expressed in terms of Goncharov's multiple polylogarithms [66J; see also |3S]. Moreover, 
relations of infinite S-sums with Xi being roots of unity have been explored in [5] . 

In this article we aim at extending these results further to important technologies that are 
currently available for harmonic sums. First we observe that S-sums can be expressed in terms 
of a variant of multiple polylogarithms involving also the parameter n; sending n to infinity 
gives back the standard generalized polylogarithms. As a consequence, S-sums can be ana- 
lytically continued up to countable many poles. This gives rise to introduce Poincare iterated 
integrals [103, 83, 45J which we call generalized harmonic polylogarithms and which extend the 
harmonic polylogarithms from [105] . It will turn out that an extended version of the Mellin 
transform applied to weighted generalized polylogarithms can be expressed in terms of S-sums. 
Conversely, the inverse Mellin transform enables one to express a certain class of S-sums in terms 
of weighted generalized polylogarithms. In addition, we extend the ideas of [1U5] to calculate the 
series expansion of generalized polylogarithms where the coefficients of the expansion are given in 
terms of S-sums. This enables one to transform and back-transform generalized polylogarithms 
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evaluated at constants in terms of infinite S-sums. In addition, we obtain algorithms to calculate 
the differentiation of S'-sums and calculate the asymptotic expansion of such sums. In general, 
we obtain complete algorithms for S-sums (1.3) with Xi E [—1, 1] and Xi ^ 0. For arising calcu- 
lations in QCD we succeeded in extending the algorithms for larger ranges of Xi by using, e.g, 
the integral representation of S'-sums. Finally, we calculate algebraic, structural and duplication 
relations for S-sums, e.g., for the alphabet X{ E {— |, i, 1, — 2, 2} to assist the compactification 
of S-sum expression. This alphabet is of special importance in case of heavy flavor corrections 
to deep-inelastic scattering at the 3-loop level. Results for related alphabets are obtained anal- 
ogously. In addition, inspired by the work in Ref. |28| for multiple zeta values we deal with new 
relations for infinite S-sums coming from the relations of finite sums, shuffle relations from the 
generalized polylogarithms and duality relations based on the argument transformations of the 
generalized polylogarithms. 

All the underlying algorithms have been implemented in J. Ablinger's Mathematica package 
HarmonicSums. For the available commands illustrated by examples we refer to Appendix 13 see 
also [B]. 

The outline of this article is as follows. In Section [2] we define the basic properties of S-sums 
with their integral representation and the closely connected generalized polylogarithms. In Sec- 
tion [3] we state important argument transformations of generalized polylogarithms that will be 
used throughout this article. In Section [4] we show how one can calculate the coefficients of the 
series expansion of generalized polylogarithms, and in Section [5] we use this result to derive an 
algorithm that transforms in both directions generalized polylogarithms evaluated at constants 
and S-sums at infinity. In Section [6] we elaborate new aspects on the analytic continuation of 
S-sums and the special case of harmonic sums. In Section [7] we define an appropriate Mellin 
transform for generalized polylogarithms and derive algorithms that calculate the Mellin trans- 
form in terms of S-sums and that perform the inverse Mellin transform in terms of generalized 
polylogarithms. In Section [8] we show how one can differentiate S-sums with the (inverse) Mellin 
transform or with the integral representation. In Section [9] we explore the algebraic relations in- 
duced by the quasi-shufne algebra, the structural relations based on the differentiation of S-sums 



and duplication relations. In Section 10 we consider relations of infinite S-sums. Finally, we 



present techniques to calculate the asymptotic expansion of S-sums in Section 11. In conclusion 



we discuss in Section |i~2| examples of generalized harmonic sums emerging in massless and massive 
3-loop calculations and their representation in the complex n-plane. In the Appendix we sum- 
marize the available commands of the HarmonicSums package which supplement the presented 
algorithms. 



2 Basic properties of 5-Sums and generalized polylogarithms 

Throughout this article let K be a field containing the rational numbers Q, let M + be the positive 
real numbers, and let N = {1, 2, . . . } be the positive integers. For a>i E N ancj^jxj E K* we define 
the class of generalized harmonic sums |90j . in short they are called S-sums, by 

n x h h x i2 ik ~ 1 x ik 
S ai ,...,a k (xi, ■ ■ ■ ,x k ;n) = ■■ ■ (2.1) 

11=1 1 12=1 «fc = l 



for a different writing see (1.3). k is called the depth and w = Y2i=o °» * s cau ed the weight of the 
S-sum S au ... Ak (xi, ...,x k ;n). 

Restricting the S-sums to the situation xi = ±1 yields to the well understood class of harmonic 
sums |19 }ll23j . For harmonic sums the limit n — > oo of (2.1 ) exists if and only if a\ ^ 1 or x\ = — 1; 
see |139j . For S-sums this result extends as follows; for a detailed proof see [6]. 

3 For a set A we define A* := A \ {0}. 
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Theorem 2.1. Let a\,a2,---a k G N and x\,x 2 ,...x k G R* for k G N. T7ie S-sum 
Sa 1 ,a 2 ,...,a k (%i,%2, ■ ■ ■ ,Xk',n) is absolutely convergent, when n —> oo, «/ and only if one of the 
following conditions holds: 

1. \x\\ < 1 A |xix 2 | < 1 A ... A \x\X2 ■ ■ ■ Xk\ < 1 

2. ai > 1 A |xi| = 1 A |x 2 | < 1 A . . . A |x 2 • • • < 1. 

7n addition the S-sum is conditional convergent (convergent but not absolutely convergent) if and 
only if 

3. a\ = 1 A xi = -1 A |x 2 | < 1 A . . . A \x 2 ■ ■ ■ x k \ < 1. 

If a S-sum is finite, i.e., the limit n — > oo exists, we define 

Sai,... ) a fc (a?l, • • • ! x k ; oo) := lim S ai a . (xi, . . . , x fc ; n) . 

An other crucial property is that for the product of two S-sums with the same upper summation 
limit the following quasi-shuffle relation holds (see |9U|): for n G N, 

S i,..,a fe (^l, ■■■,x k ;n) S 6ll ...,6 ! (yi, ...,yr,n) = 

n j 

^ ^S a2i ... jafc (x 2 , ...,x k ;i) S 6l) ... i6i (yi, ■■■,yui) 



i u 
i=i 



r b\ 



+ ^2 -^[\Sai,...,a k (xi, ...,x k ,i) S b2i ... ife; (y 2 , ■■■,yr,i) 

i=l 

~ X] S a 2 ,-,a fc (^2, • • • , Xfe, ») 6,(lfti, • • ■ , VTi • ( 2 - 2 ) 



i=l 



Since the right hand side consists of S'-sums where at least one component of the product 
has smaller depth, this formula applied iteratively leads to a linear combination of S'-sums with 
integer coefficients. The underlying quasi-shuffle algebra induced by this product property will 
be exploited further in Section |9| 

Finally, we will exploit heavily the fact that each S-sum can be written as an indefinite nested 
integral. Namely, for Si(6;n) with b G M* we have 



x n — 1 

Si(6;n) = / — -dxi 



o xi - 1 



and by induction on m > 1 we obtain 



S m (6;n)= / — / / — / — -dx 1 dx2 ■ ■ ■ dx m -idx m . 

Jo x m Jq x m -\ Jq X2 Jo x± — l 

Now consider, e.g., the sum S 2 i{a, b; n) with a, b G K*. Due to the integral representation for 
S-sums with depth 1 we obtain 



lb A 

^ — ^ qJ 

S 2 ,i(a,6;n) = 2^^Si(b;j) 



Aa^ r b xi -1 
j^P J x-1 X 

1 n j 
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(S2(a x; n) — ^{a; k)) dx 



o x - 1 

1 /• <BB 1 r y z n -I 
x-1 J a y J z-1 
ab ^ r x 1 r y Z n — 1 



x-a J a y Jo z-1 



dzdydx 
dzdydx. 



Looking at this example one arrives straightforwardly at the following theorem; the proof follows 
by induction on k and the integral representations of S-sums with depth k = 1. 

Theorem 2.2. Let rrii e N, bi £ R* and define 



Irn\ ,ni2, ...,mfc (^1) ^*2> • • • j %) — 

rb v -b k dx rn k rfx™^ 1 f X l dx\ f X t d X \ 

Jo <Wo xf^ 1 "Jo ~4Jo x\-b 1 ---b k . 1 

Jbvbk-! Xk-i 1 Jo xT-l 1 ' 1 J° X k-l Jo 4-l" 6 l---^-2 

r*l-i dx™*~ 2 f*C- 2 ~ 2 dx™*^ 1 [4-2 dx\_ 2 f<-2 dx\_ 2 

Jb 1 -b k _ 2 ^-2 2 JO ^0 X *-2 A 4-2- 6 l---^-3 



dx™ 3 f x 3 dx™ 3 ' 1 /* X 3 ^2 ,xg ^1 



bi& 2 6 3 ^s 13 ./o x™ 3 1 7o a?| JO 4- fe i & 2 



™2 / m 2 -l 
bl6 2 X 2 •/ X 2 JO 



/ 

•76 



a dx™ 1 p dx™ 1 " 1 dx? ^ (x 1 ) 2 - 1 ! 

(XX i . 



T m i / mi-1 / „2 / ™1 i 1 

bi x l JO X-y Jo x l Jo x l — 1 

T/ien /or a// n G N, 

Smi,m 2 ,...,m fc (^lj &2> ■ ■ • , ^fc! n ) = ^mi,m 2 ,...,m fc &2> • • ■ ,&fe! n )- 

Note that sending n — >■ oo we obtain back multiple polylogarithms stated, e.g., in [SSI E3 EO] ■ 

For later constructions based on analytic reasoning (see also Section [6]), we will use the fact 
that one can transform S-sums to a particular form of integral representations, the so-called gen- 
eralized harmonic polylogarithms using the inverse Mellin transform; see Section [7J Conversely, 
with the Mellin transform one can express a weighted generalized harmonic polylogarithm in 
terms of S'-sums. In other words, the S'-sums are strongly twisted with generalized harmonic 
polylogarithms which can be defined as follows. 
For a£i and 

a, if a > 
oo, otherwise 



we introduce the auxiliary function f a : (0, q) \— > M by 



fa(x) 



1 

| a. | — sign(a) x ' 



if o = 
otherwise. 



(3 



Then the generalized harmonic poly logarithms (in short generalized polylogarithm) H of the word 
mi, . . . , m w (resp. vector m = (mi, . . . , m w )) with mi S 1 is defined as follows. For x G (0, <?) 
with g := min m,, 

nii>0 

H(x) = 1, 

i(logx)^, if (m 1 ,...,m w ) = {0,...,0) 



H M / if (mi,... 

ti mi ,m 2 ,...,m w [x) - | Jq 31 f mx (y)Hm 2 ,...,m w (y) dy, otherwise. 



Generalized polylogarithms are also dealt with in Refs. [391 113 EH El ED] as hyperlogarithms. 
The length w of the vector m is called the weight of the generalized polylogarithm H m (x) . 

Restricting Hm 1 ,m 2 ,...,m w (x) to rrii G {—1,0,1} gives the class of harmonic polyloga- 
rithms |105j . Typical examples are 

f x 1 

Hi(x) = / dxi = - log(l - x) 

Jo 1-xi 

H_ 2 (x) = f ^—d Xl = log(2 + x) - log(2) 
7o 2 + xi 

^A-l^) = / / — / TT — dx 3 dx 2 dxi. 



x\ Jo x 2 Jo J + x 3 



A generalized polylogarithm H m (x) = H mi ,...,m w (x) with q := min mi> omj is an analytic 
function for x G (0, g). For the limits x — )■ and x —> q the following holds. From the definition 
we have that H m (0) = if and only if m ^ W . In addition, H m (q) is finite if and only if one of 
the following cases holds 

• mi 7^ q 

• w > 2, mi = 1 and m„ = for all v with 1 < v < w. 

We define H m (0) := lim x ^ + H m (x) and H m (l) := lim ;c _ s>1 - H m (x) if the limits exist. 
Note that for generalized polylogarithms of the form H mii ... jmfej i j o...,o(^) with c := min mi> omj 
(and c > 1) we can extend the definition range from x G (0,1) to x G (0,c). These generalized 
polylogarithms are analytic functions for x G (0, c). The limit x — > c exists if and only if mi ^ c. 

Summarizing, a generalized polylogarithm H mi ,m 2 ,...,m w ( x ) with q : = min mi> omj is finite at 
c£l for c > if one of the following cases hold: 

1. c < q 

2. c = q and mi / q 

3. (mi, . . . , m^,) = (mi, . . . , m^, 1, 0, . . . , 0) with := min mi> o m« and c < q' . 

Throughout this article, whenever we state that H mijjn2) ... imiu (c) is finite or well defined, we assume 
that one of the three conditions from above holds. 

For the derivatives we have for all x G (0, min m .>o m,j) that 

~, Hra(^) — /mi (^)Hm2,m3,...,m„,(^) ■ 

dx ~ 

Besides that, the product of two generalized polylogarithms of the same argument can be ex- 
pressed using the formula 

H p (x)H ? (x)= R r(x); (2.3) 

r=pLLlq 

here the shuffle operation pmq denotes all merges of p and q in which the relative orders of 
the elements of p and q are preserved. Thus, similarly to iS-sums, one can write any product 
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of generalized polylogarithms by a linear combination of generalized polylogarithms over the 
integers. 

As utilized in |105j for harmonic polylogarithms, we use this shuffle product to remove trailing 
or leading constants in the vector of the arising generalized polylogarithms. E.g., for a£l the 
shuffle product gives 

Ho(in) H)7i lj .„ ( j7i 11l \X) — Htt,mi,...,m m (jE) ~\~^mi,a,m2,...,m w {,x) 

+ H mi) m 2) a,m3,. ..,rn m {%) + " ' ' + H mij ... )7nt(J)a (ic) , 

or equivalently 

^■mi,...,m w ,a{%) — H a (x) H mii ... : rn m H a 7Tll j jJ71ll) (x) 

Hmi,a,m2,...,m t |,(^) ' ' ^mi,,..,a,m w (p') • 

If m w is a as well, we can move the last term to the left and can divide by two. This leads to 

H-mi,...,a,a(^') — ^ (-^-a^) ^wii, ...,m w -\,aifl) Ha,ini,...,m„_i,a(^) 

Hmi,a,m2,...,ra al ^i,a(2') ' H-mi ,---i a i m iu-i (*^)) " 

Applying this tactic iteratively leads to a polynomial representation in H a {x) with coefficients 
in terms of generalized polylogarithms with no trailing a's. In particular, if a generalized poly- 
logarithm is analytic for (0, q) then also the occurring generalized polylogarithms of the linear 
combination are analytic for (0, q). A typical example for a = is, e.g., 

H 2) - 3> o,o(x) = ^Ho(x) 2 H 2 ,- 3 (x) - H (x)Ho,2,-3(x) - H (x)H 2 ,o,- 3 (z) 

+Ho,Q,2,-3(aj) + H ,2,o,-3(aj) + H 2 ,o,o ,-z(x). 

Note that there is one extra complication: given a generalized polylogarithm in the form 
H-m 1 ,...,m k ,i,o...,o(x) with 1 < q := min mi> om,j then we cannot use the strategy to extract trailing 
zeroes if x > 1 since this would lead to infinities: The shuffle algebra would suggest to rewrite 
£[2,1,0(20 being analytic for x G (0,2) as Ho(x)H 2j i(x) — Ho, 2 ,i(a;) — H 2) o,i(x) , however for x > 1 
for instance H 2j o 5 i(x) is not defined. In later considerations these difficulties have to be handled 
separately. 

Completely analogously, leading constants can be removed in generalized polylogarithms. In 



addition, using the quasi shuffle algebra (2.2) it is possible to extract leading (or trailing) indices 



of a S'-sum. For example we have for a G N and b G K* 

S a ,2,i [b, 5, §;n) =S a (6;n)S 2 ,i (§, |;n) + S 2 +a,i (§, §;n) + S 2 ,i+ a (§, 
- S 2 ,a,i (5, b, |;n) - S 2) i ia (A, |, 6; n) . 

For 5-sums S ai) ... )afe (&i, . . . , 6^; n) with Oj 6 N and 6j G [—1, 1] \ {0} this can be used to get a 
polynomial representation in Si (re) with coefficients in terms of convergent S*-sums, like, e.g., 

S ljl>2 (1, 1, |; re) = ±Si(n) 2 S 2 (|;n) + Si(re) (S 3 (§; n) - S 2> i (A, l;n)) 

- S 3 ,i (|, 1; re) + S 2j i,i (§, 1, 1; n) . 

For general S'-sums this strategy usually fails to extract the divergency, e.g., in 

Si, 2 ,i (5, g, g;n) =Si(5;re)S 2 ,i (§, + S 2 , 2 (|, |;n) + S 3) i (|, |;n) 

- S 2 ,i,i (5, |,5;re) - S 2 ,i,i (|,5, |;re) , 

the sums Si(5;n),S 3i i (§, |;re) and S2,i,i (5 ,5, |;n) are divergent. 



8 



3 Identities between Multiple Polylogarithms of Related Argu- 



In the following we want to look at several special transforms of the argument of generalized 
polylogarithms which will be useful for later considerations. In general, we start with a defined 



x € M and derive representations in terms of generalized polylogarithms in x. Since such a 
transformation holds not only for a point but in a region around x, all the transformations 

have the following additional property: if the input expression f(x) is analytic in a certain 
region around x, also the arising generalized polylogarithms of the transformed representation 
are analytic in this region. 



3.1 x + b — > x 

Lemma 3.1. Letmi 6 R\(0, 1) andx > such that the generalized polylogarithmH mit ... tmi (x + 1) 
is defined. If (mi, . . . , m/) / (1, 0, . . . , 0), 

Hmi,...,m; [x + 1) =H mi) ... )mi (l) + H mi _i(x) Hm2,...,mi (1) + 

+ • • • + H mi _i ) ... )mj _ 1 _i(ar) H m( (l) + H mi _i v .. ;m( _i(x) . 

7/(mi,...,m,) = (l,0,...,0), 

Hmi,..., mi (x + 1) =Hi A ..., (l) - H _i,... _i(x) . 

Proof. We proceed by induction on I. For I = 1 and mi > 1 we have: 



ments 



generalized poly logarithm f(x) = H. 




with particular chosen integers a, b, c, d and 




For I = 1 and mi < we have: 




For I = 1 and mi = we have: 




Suppose the theorem holds for / — 1. If mi / Owe get 



H. 
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"l"H jrl2 _i(x) H m3ij ,^(1) + H m2 _x iTra3 _i(x) H-rn4,...,mi (1) ~i~ 

+H rrV2 _i j ... jmi _ 1 _i(a;) H mj (l) + H m2 _i v .. ;m( _i(x)^ dn 

— Hmi,...,mj(l) "I" H mi _i(x) Hjj l2 ,...,m; (1) 

+H mi _i !m2 _i(x) H m3v .. ;m( (l) + • • • + 
+H mi _i ) ... jmi _ 1 _i(a;) H mj (l) + H mi _i v .. )m( _i(x) . 

If mi = we get 

H , m2 ,..., m( (x + 1) = H , m2 ,..„ m ,(l)+ T +1 Hm2 -- m ' ( " ) d7x 

- H , m2 ,..., m; (i) + y o — — dn 



- H , m2! ..., mi (l) + ^ 



m2,...,mj 



(1) + 



+Hm2-l( 2 ') -^"13, ••■, m ! (-0 "I - ^ni2- l,m3— lO^) H m4i ... im; (1) + 

H h H m2 _i ) ... jmj _ 1 _i(a;) H m; (l) + H m2 _i v .. jm( _i(x)^du 

= Ho, m2 ,..., mj (l) + H_i(ar) H m2j ... jmi (l) 
+H_i 

,m 2 -l( x ) H m3) ... )mi (l) + • • • + 

+H_i,..., mi _ 1 _i(x) H m; (l) + H_i,... )ro ,_i(x) . 



□ 



The proofs of the following three lemmas are similar to the proof of the previous one, hence 
we will omit them here. 

Lemma 3.2. Let b > 0, m\ G R\(0,6] ; m, G R\(0,6) fori G {2,3, . . . 1} such that (rrij, ...,mi)^= 
(1,0, ...,0) for all j G {1,.../}. IfH mu ... jm (x + b) is defined for x > 0, 

Hj7i lj ... j m i (aj + b) = H mijm2i _ m; (6) + H mi _fe(x) H m2) _ irrt; (6) 

~t~Hmi— 6,m 2 — &(-^) H m3i . . m; (o) + • • • + 

~l~Hmi— b,...,mj_i— bO^) H m; (6) + Hj^—^ !m( _;,(x) . 



Lemma 3.3. Xei (mi, . . . ,m/) = (1,0, ... ,0) and let x > such that H mi) ... ;m; (x + 5) is defined. 
For 1 > 6 > and 1 - b > x > 0, 

Hmi,...,m; (x + 6) = H mijJ7l2i _ i?Tl! (6) + H mi _f,(x) H m2 _ irn; (&) + H mi _^ m2 _^(x) H m3i ,..^(6) 

(x). 

For x > ana 1 6 > 1, 

Hmi,...,mj(^ "I" fr) = Hmi,m.2vi m ((^) H mi _j(x) H m2) __ jmi (6) H?rti — 6,m 2 — b(^) H m3i ... )?Tt; (6) 
... H mi _fe j ... 5 m; _i— b(^) H m; (6) H mi _b ) ... ) m i — ft(x) . 



Lemma 3.4. Zei 6 > 0, mi G R \ (0, 6], m,eM \ (0, 6) for i G {2, 3, ... k} and (m k+1 , . . . ,m;) 
(1,0,..., 0). Xei x > snc/i that H mi; ... im; (x + 6) is defined. For 1 — 6 > x > and 1 > 6 > 0, 

Hmi,...,mi{% + 6) = H mij?Tt2j ... im; (&) + H mi _ft(x) H m2v .. i?ni (6) + • • • + 
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~i~ ^-mi—b,...,mi—i— H m; (6) -|- H mi _b v .. ) m(— ft(x) • 



For 6 > 1, 
H 



mi,...,mj(^ "I - 6) — H mi l m,2,...,m; (&) "I - H-mi — b(x) H m2 ,...,m; (^) "I - ' ' ' 

+ H mi _6 r .. jmfc _6(x) H mfc+lj „. : mi \y) ~ Hm 1 -b,...,m k+ i-b{%) H mfe+2r ,. jTOz (6) + 

H 1- H mi _6 ) ... )m j_ 1 _6(a;) H m; (6) — H mi _ fcj ... iTOj _&(x) 



3.2 b — x —7- a; 

In this Subsection we assume that b > and we consider indices £ R \ (0, b). The following 
method transforms a defined generalized polylogarithm ff roii ... lTO „(i — x) for some x G R in 
terms of generalized polylogarithm with argument x and evaluations at constants. Proceeding 
recursively on the weight w of the generalized polylogarithm we have 

H mi (& - x) = H mi (6) - R b _ mi (x) for mi/ b, 
R b (b -x) = H (x) - Ho(6) . 

Now let us look at higher weights io > 1. We consider H mit m 2 ,...,rn w {b — x) with m; G 1\ (0, b) 
and suppose that we can already apply the transformation for generalized polylogarithms of 
weight < w. If mi = b, we can remove leading b's (see Section [2]) and end up with generalized 
polylogarithms without leading b's and powers of Hb(b — x) . For the powers of Hfe(6 — x) we can 



use (3.1); therefore, only the cases in which the first index mi / b are to be considered. For 
6/0: 

H/t \ _ tt / t \ _ ^ Hm2,...,m m (fo ~ Q 11 
mi ,r7i2,...,mu, XJ — -H-mi ,m2,...,m» / | , zj „: /" _ , lW 



-mi + b\ — sign (—mi + b) t 

Since we know the transform for weights < w, we can apply it to H m2j ... jm , m (6 — t). Linearizing 
the expression with the shuffle product using the definition of the generalized polylogarithm we 
finally obtain the required weight w identity. With the steps above one obtains, e.g., 

H 3 ,_ 2 ,-i(2 - x) = -H_i(x)H_ 2) _i(2) + H_i(2)H_ M (x) - H_i i4)3 (x) + H 3 ,_ 2 ,_i(2). 

Remark 3.5. Analyzing the underlying algorithm we see that in the transform of 
(b - x) we find only one generalized polylogarithm of weight w and argument x, 
namely H_( mi _ b ) j _( m2 _ fo ) ) ...._( m „_b)(x). All other polylogarithms at argument x that appear have 
lower weights. In general the indices of the generalized polylogarithms at argument x that appear 
in the transformed expression are in the set { — (mi — b), — (m 2 — b), .... — (m w — &)}. Note that 
if we look at the reverse transform i.e., x — > b — x we can derive a similar property. 

3.3 ~, — ~ y x 

l+x 

Next, we present an algorithm that performs the transformation — > x, like, e.g., 

H-2,i (r^j ^-iW (-H-s(x) + H_i(x)) + H_ 3 ,-i(x) - H_ 3 ,o(x) + H_ 2 ,i(l) 
-H_!_i(x) + H_i i0 (s). 
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Here we consider indices with m; € 1\ (0, 1). Proceeding recursively on the weight w of the 
generalized poly logarithm we have for < j=§ < 1, yi < — 1, — 1 < y% < and y% > 1 : 

1 — x s 



H_i 

H 
Hi 



1 + x 
1 — X 
1 + x 
1 — X 
1 + x 
1 — X 
1 + x 
1 — X 
1 + x 
1 — X 
1 + x 



-Ho(-yi) - H_i(x) + H (l - yi) + Hi^(x) 
H_i(l) - H_ x (x) 

-Ho(-ite) - H_i(x) + H (l - y 2 ) - Hi- V2 (x) 
-Hi(x) + H_i(s) 
-H_i(l) - Ho(x) + H_i(x) 



H (y 3 ) + H_i(x) - H (y 3 - 1) - H^W 

l+»3 



(3.2) 
(3.3) 



Now let us look at higher weights w > 1. We consider H mi)m2) ... )TOlu ^ j+fj with m, 6K \ (0, 1) and 
suppose that we can already apply the transformation for generalized polylogarithms of weight 
< w. If mi = 1, we can remove leading ones and end up with generalized polylog arith ms without 
leading ones and powers of Hi(;[=f). For the powers of Hi I yrf ) we can use rt3.3h; therefore, 
only the cases in which the first index mi ^ 1 are to be considered. For y\ < U, y\ ^ —1 and 
2/2 > 1 w e get: 



H 



yi,m,2,...,m v 



H— l,m2,...,m„ 
ho,m2,...,m m 



H 



j/2,m 2 ,...,m„ 



1 — X 

1 + x 



1 — x 
1 + x 
1 — x 
1 + x 



1 — x 
1 + x 



1 



1 + t 



-H 



ra!,...,m, 



1 -t 



dt 



+ 



H_ 



1-t 



-H 



dt 



l,m,2,...,m, 



,(1) 



1 - t 



H r 



dt 



Ho,m 2 ,...,m TO (l) 

1 H 

r-tl r . 



1-t 



i - r m2 '-' m ™ v i + 1 



-H r 



dt 



1 + r - 2 ,-,m 1 „ u + t 



1 - t 



dt 



H-y2,m2,---,m w (l) + 



-H r , 



o 1 + * 



i^ ^ ™» u + t 

1+2/2 



rH m2i m , 
1-t 



1-t 
1 + t 



dt 



dt. 



(3.4) 



(3.5) 



1+x 



1+x 



Note that the transformation (3.4) can be only applied if (mj,- • • ,m w ) is not the zero vector. 
However, by definition the arising special case 'Hq > ....()( tt^\ reduces to Hof rr?l which we handled 
in fliT2 ). 

Since we know the transform for weights < w, we can apply it to H m2) ... jTOtu ( £pf J • Linearizing 
the arising products of generalized polylogarithms with the shuffle product and using the definition 
of the generalized polylogarithms we obtain the required weight w identity. 



3.4 kx — >■ x 



Lemma 3.6. If nn G R, m\ ^ 0,x € K + and k £ M + such that the generalized poly logarithm 
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H mi ... m Ak ■ x) is defined then we have 



H-mi,...,mi(k ■ x) = Hmi ^ mi(x) . (3.6) 



fe k 



Proof. We proceed by induction on I. For / = 1 we have: 

rkx i t' x k 

H mi (fc-x) = / | . . — —du = I | . . — — — du = Hmi(aj) . 

Jo |mi| — sign(mi) u Jo l m il ~~ sign(mi) few, * 

Suppose the lemma holds for I. For mi / we get 

XT /, \ f x H m2 . mi+1 (u) H m2v .. jm ; (/cn) 

H mi ,... )TO , +1 (fe • x) = / -J . . — —du= / . . — — — kdu 

Jo \mi\ — sign(mi) u J \m\\ — sign (mi j ku 

j-X Hmj m i + 1 (u) 
Jo J -^ Li — sign (mi ) u k ' ' k 

For mi = we get 

tt /, \ f X Hm2 ,— ,mi+\ ( u ) j f x H m2 ... m , (feu) 
Ho,ma,...,mj + i(« = / ~ du = j — kdu 



l>X Hmj m i+i (u) 

/ fc — k du = Ho rn2 ™!+l(x). 



□ 



The proof of the following lemma is straightforward (compare |105j . where it is stated for harmonic 
polylogarithms) . 

Lemma 3.7. Ifnrii G M,mi ^ andx £ M - such that the generalized polylogarithmH mit ,.. tmi (—x) 
is defined then we have 

H mij ... jmi ( x) = ( 1) H_ mi) ... i _ JTli (x) , (3-7) 
feere fe denotes the number of m, which equal zero. 

3.5 

Here we assume that < x < 1. First, we consider only indices with m-j < 0. Proceeding 
recursively on the weight w of the generalized polylogarithm we have for y < that 

H y f -) = H i (x) - H (x) - H (-y) and H ( -) = -H (x) . 



Now let us look at higher weights w > 1. We consider H mi>m2j ... im f-J with m; < and suppose 
that we can already apply the transformation for generalized polylogarithms of weight < w. For 
mi < we get (compare [105J ) : 

Hmi,m2 ""(i) = R ^™ 2 ,...,mM) + ^ t 2(_ mi + ^.-."^ ( j ) d * 

= H mim2j ... )TOto (1) + J — H m2i .. mu ( — I c?t 



1 i /l 

Hm.Q nx... I — I dt 



mi 
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/1\ f 1 1 fl\ 

Ho,m 2 ,...,m M y — J = Ho,m 2 ,...,mi„(l) + y t 2 H /t)^ mj '" ,,m " ^ 

= Ho im2) ... )TOw (l) + y -H m2 ,..., mM T-^ dt. 

Since we know the transform for weights < to, we can apply it to H-m 2 ,...,m w (j) and finally we 
obtain the required weight w identity by using the definition of the generalized polylogarithms. 

We can extend this transformation by the following observations. An index m > in the 
index set leads to a branch point at m and a branch cut discontinuity in the complex plane for 
x G (m, oo). This corresponds to the branch point at x = m and the branch cut discontinuity 
in the complex plane for x G (m, oo) of log(m — x) = log(m) — H m (x) . However, the analytic 
properties of the logarithm are well known and we can set for < x < 1 for instance 

H m (-) = K ± (x) + R Q (m) + R (x) -iir, (3.8) 

\X J m 

by approaching the argument - form the lower half complex plane. We can now consider an 
alphabet with at most one letter m > 1 (Note that we could as well assume m > 0. However, for 
simplicity we restrict to m > 1). For such an alphabet the strategy is as follows: if a generalized 
polylogarithm has leading m's, we remove them and end up with generalized polylogarithms 
without leading m's and powers of H m (-) . To deal with the generalized polylogarithms without 
leading m, we can slightly modify the previous part of this Section and for the powers of H m Q) 



we can use (3.8). Following this tactic, we arrive, e.g., at the identity 

H-2,1 (-) =-t7r (-H_i(l) + H_i(x)-H (x)) +H_ 2)1 (l)-H_i (l) + H_i (x) 

\X J V 2 2 / 2' u 2' u 

- H_i + H_ ± x {x) - H ,o(x) + H ,i(l) - H ,i(x). 

2' 2> 

4 Power Series Expansion of Multiple Polylogarithms 

An algorithm to calculate series expansions for harmonic polylogarithms Hm 1 ,m 2 ,...,m w (x) with 
mi G {—1,0,1} has been presented in |105j . Subsequently, we generalize these steps for gener- 
alized polylogarithms, i.e., m, G M. In general, H TOl)m2j ... jmTO (a;) does not have a Taylor series 
expansion if trailing zeroes occur, i.e., the expansion consists of a log(x)-part and a part free of 
log(x). The separation into these two parts can be easily accomplished by removing trailing zeroes 
as described in Section [2j In other words, one gets a polynomial expression in Hq(x) = log(x) 
with coefficients in terms of generalized polylogarithms that have no trailing 0's. For each of 
these coefficients we are now in the position to calculate a Taylor series expansion as follows. 
For polylogarithms of weight 1 with j G R and x G (0, \j\) we have 



Rj(x) 



Let m = (mi, . . . , m w ) be a non-zero vector and q = min m -^o Assuming that for x G (0, q) 



oo 

H m (x) = ^^r s "( b;i) 



i=i 



*For a definition of single- valued polylogarithms see [38l [5T1 H9l |44] . 
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one can write the relations for j > and x G (0, min(j, q)) 

CO j j 

H , m (x) = E^+T S "( b;i ) 



oo 



IT 

i=l J 

oo ,■ 
X 



E 777 s ct,«(crj, 6; t) - £ ^S„(6; i) 



=1 U i=l * 



i=l H J; 



00 



i=l v J; i=l 
Proof. We just prove the case that j > 0. The other cases follow analogously. 

ttjA x ) = / — 2^-^r S n(b;i)dy 
Jo J y i=1 » 

„.£ -. OO / \ fc OO j 

= 1 / E(-) E^kV s ^ + 1 ^ 



fe=0 w y i=0 



/-x °" * / \ »-« fc+1 fe+1 

- U gg(f) w^*™ 



00 ,■ 



= E^v^' 6 ''- 1 )- 



j=i 



□ 



Separating the log(x) = Hq(x) part from the log-free part and applying the formulas from 
above we get, e.g., 

it , , y 2-*(-x)* " 2"*(-x)* 



Summarizing, given H mijm2) ... jrriu ,(ar) with m w 7^ and q = min m ^o \ m i\i one can calculate 
a Taylor series expansion which holds for x G [0, <?) and where the coefficients of the Taylor 
expansion are given explicitly in terms of 5-sums. 

Inverse construction. In addition, since the formulas from above can be reversed, also the 
other direction is possible: given a sum of the form YlTi^Y^r 1 with & G N U {0} and c G M*, 
one can compute a linear combination of (possibly weighted) generalized polylogarithms h(x) 
such that for x G (0, j) with j > small enough YlT( cx T^v^ = h ( x )- 
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E.g., for i G (0, |) we have 



E 

i=i 



2VSi,i (3,i;i) 



H n i (x) + H n i i (x) + Hi i (x) + Hi i i (x). 

u ' u ' ft U !R>S ' 9' u 'fi J1BH 



Asymptotic behavior. Moreover, one can determine the asymptotic behavior of a generalized 

' ( 1 M M .11 VI . >.- I (HI II vv, 1 

y 



polylogarithms H m (x): Define y := ~. Using Section 3.5 on H m (^-J = H m (x) we can rewrite 

H m (x) in terms of generalized polylogarithms at argument y together with some constants. Now 
we obtain the power series expansion of the generalized polylogarithms at argument y around 
easily using the previous part of this Section. Since sending x to infinity corresponds to sending 
y to zero, we get the asymptotic behavior of H m (x) . 
E.g., performing these steps yields 



H_ 2 , (x) =H_ 2 ,o(l)+H_i i0 (l) 
=~ H (x) 2 - H (x) 




-2,0 



1) +H 



)(!)• 



Values of Multiple Polylogarithms Expressed by 5-Sums at 
Infinity 



Consider the generalized polylogarithm H mi 



m2,...,m v 



(x) with q := min mi> o ntj and c£R, c > 



such that H 



mi ,rri2,...,m w 



(c) is finite. To be more precise, one of the following cases holds: 



1. c < q 



2. c = q and mi ^ q 

3. (mi, . . .,m w ) = (mi, . . . , m k , 1,0, ... ,0) with q' := min mi>0 mj and c < q' . 

If the generalized polylogarithm H m (c) has trailing zeroes, we first extract them. If < c < 1, 
we end up in an expression of generalized polylogarithms without trailing zeroes and powers of 
Ho(c). Otherwise, if c > 1, we cannot handle polylogarithms of the form H miim2j ... jmfej i o,...,o(c) 
(see Section [2]) and keep them untouched. Note that by the shuffle product no extra indices are 
introduced and the first non-zero index is again m\. Thus the arising polylogarithms are again 
finite at the evaluation point c. 

We can use the following lemma (iteratively) to rewrite powers of Ho(c) . 

Lemma 5.1. Let c > 0. We have 



< c < 1 
c = 1 
1< c < 2 

2 < c 



H (c) 
H (c) 
H (c) 



-Sifl 



-Sifl 



c; oo 



c; oo 



H (c) = Ho(2)-H ^- 



To rewrite generalized polylogarithms of the form H mi . m2v .. jjnm (c) where m w ^ and 
mhim.^o \rrii\ > c, we use the power series expansion of H mi ,m 2 ,...,m, w (x) and send x — > c. E.g., for 
one of the arising sums we obtain 



i=l 



Y"] °-^S n (b; i) ->■ S tt) „(c<r, b; oo) 
L — ' i a 
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What remains open are the situations of the type H mi , m2 ,..., mfc ,i,o,...,o(c) with c > 1 and 
polylogarithms H mi , m2 ,..., mt() (c) with m w / 0, max mj <omj > -c, m.m mi>0 mi > c and mi / c. 
As shortcut, on© c&jH rewrite H 7ni5?Ti 2 ) ... 5 m ltJ 

(c) with (mi,ro2, . • • , m^) = (1, 0, . . . , 0) and c > 1 

by using the following lemma iteratively. 

Lemma 5.2. For c > 1 and (mi, m 2 , . . . , m^) = (1, 0, . . . , 0) we /iave 

H mi , m2 ,..., mffi (e) = (-l) w+1 S m (l; oo) - H , mi _i,..., mu ,_i(c - 1) . 
To rewrite generalized polylogarithms of the form H mi , m2 ,..., mro (c) where mi G E \ (0, c], 



m w 7^ 0, min mi> omj > c and s := max m .<omj > — c we use Lemma 3.2 and afterwards send 
x — > c + s, i.e., we get 

H mii ... )?Ttw (c) =H m . 1 , ni2 ,..., mtu (— s) + H nj , 1 _|_ s (c + s) H rn2) ,.. )niTO (— s) 

+ H mi _|_ SjJrl2 _|_ s (c + s) ^m 3 ,...,ra w { — s) + h 

(c + s) H mi (-s) + H mi+S ,..., mm+S (c + s) . 
Note that — s < c and c + s < c. To rewrite generalized polylogarithms of the form 



Hmi,m 2 ,...,m fe ,i,o,...,o(c) where c > 1 we can proceed similar to the previous case: we use Lemma 3.4 
to rewrite H mi)m2j ... jmm (x + 1) and afterwards send x — > c — 1. Observe that both cases reduce 
the problem to a simpler situation: We can either handle all arising polylogarithms by the recipes 
given above. If not, observe that the remaining polylogarithms are of the form H ni , n , (ii) with 
q' := min ni ^o \ n i\ and u < c. Note that the evaluation point gets smaller. In addition, the index 
vector (ni, . . . ,«;) is given by a subset of (mi, . . . , m w ) up to the possible addition of s < where 
s is the maximum of the arising negative indices of (mi, . . . ,m w ). Thus applying this reduction 
iteratively (to the polylogarithms that cannot be handled by the cases from above) decreases step 
by step the positive evaluation point and in the same time decreases the negative indices and 
thus enlarges the absolute value of s superlinearly. Consequently, from a certain point on q' gets 
larger than q and the difference between these values increases in each iteration. This guarantees 
that eventually u < q' , i.e., the arising polylogarithms can be transformed to 5-sums by using 
the series expansion representation and the evaluation at x = u as described above. 

Summarizing, we can rewrite all finite generalized polylogarithms H mijm2j ... iJnro (c) in terms of 
5-sums at infinity, like, e.g., 

H 4 ,i, (3) = - S_ 3 (oo) + 2S 3 (oo) + S_ 2 (oo) Si(|; oo) - S 2 (oo) Si(§; oo) 

- S 3 (-i; oo) + S_i(oo) (S_ 2 (oo) - Si,i(i, -2; oo)) 

- Si >2 (3,4;oo) + Si )2 (g, -3;oo) + Si j2 (|, -1; oo) 

+ S 2 ,i(-1, l;oo) - S 2j i(j,4;oo) - Si ) i il (i, -2, ^;oo) . 

Note that this process can be reversed: given a finite /S-sum at infinity we can rewrite it, e.g., in 
terms of generalized polylogarithms at one: 

S 2 ,3,i( — h I' 2; oo) = — 1X0,-2,0,0,-6,-3(1) + Ho -2,0,0,0,-3(1) + Ho,0,0,0,-6,-3(l) 
— Ho,o,0,0,0,-3(1) • 

For a complete algorithm for the reverse direction we refer to [6l Algorithm 1]. 

We note that if we consider the generalized polylogarithm H miiJn2) ... jmiu (x) and c£l, c < 



such that H mi , m2 ,..., mm (c) is finite then we can use Lemma 3.7 to transform Hm 1 ,m 2 ,...,m w (x) to an 
expression in terms of generalized polylogarithms at argument — c. Hence we can again use the 
strategy presented above to rewrite this generalized polylogarithm in terms of 5-sums at infinity. 
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6 Analytic Continuation of 5-sums 



It has been remarked in |137j that Carlson's theorem \A2\ Q~22j^] can be used to perform the 
analytic continuation of harmonic sums: the evaluation of the harmonic sums at integer points 
determines uniquely (among a certain class of analytic functions with exponential growth) the 
analytic continuation in the complex plane. We will elaborate this property in further details for 
harmonic sums and more generally for S'-sums. 

A function / : M — > C for a region M C C is of Carlson- type (C-type) if 

• it is analytic in the right half plane, 

• it is of exponential type in the half plane, i.e., there exist c, C G R + such that for all z G C 
with Im(z) > Owe have that 

\m\ <Ce c W 

• and there exist c, C G M + with c < it such that 

\f(iy)\ <Ce c \y\. 



For such functions the following theorem holds; for a detailed proof (together with even more 
general statements) we refer to (33J Sec. 9.2]. 

Theorem 6.1 (Carlson). Let h(z) be a function of C-type. If h(n) = for all n G N then h = 0. 

Throughout the remaining article we rely on the following consequence: Let f{z) and g{z) be 
functions of C-type. If f(n) = g{n) for all n G N, then / = g in the domain where both functions 
are defined. 

The following simple closure properties will be used later. The proofs are immediate. 

Lemma 6.2. Let f(z) and g(z) be functions of C-type and let a G M + . Then f{z) + g{z), 
a z = e 2ln(a ) and a z f(z) are of C-type. 



Proposition 6.3. For rrii £ N, 6; G 



and z G 



consider the integral f(z) 



mi ,77i2,...,mfc 



(&i, 62, ... , bk', z) from Theorem 2.2 Then f(z) is of C-type. 



Proof. We note that f(z) is analytic for Re(z) > 0. Next we show that f(z) is of exponential 
type. For 1 > \z\ > with Re (z) > we can bound |/(^)| by a constant. For \z\ > 1 and 
d := max(&i, 6162, • • • , ^1^2 • • • 6fe) > we have 



\m\ 



bi-b k dx™ k 



< 



bi-b k dx™ k 



m k I 



dx\ 



o x\ - h ■ ■ ■ b k -i 



■1 dx™ 



• 2 2 

1 X 



-dx 



dx\ 


f x ' dx™ 1 




x z — 1 


b\ 


Jb 1 ki 1 


Jo 


x — 1 



dx. 



(6.1) 



Now we show that there is a c G M + such that 
for a G M + , b G R. If d < 1, we have 



x-l 



< c(d+l)\ z \ for < x < d and all z = a + bi 



1 



'The original theorem needs to be extended to be even applicable in case of the harmonic sums. 
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Finally, let d > 1. Since lim,^! 



x z -l 



x-l 



< (d + l)' z ', there are e > and c > such that 



x z -l 



x-l 



< c(d + l)' 2 ' for 1 — e<x<l + e. It remains to show that 



x G (0, 1 -e) U (1 + e,d) : 



< c(d + 1)I 2 I for 



1 



<^<i(d|l)N. 



In summary, we can choose c G M + (independent of the choice of z) with 
< x < d and using (6.1 ) it follows that 

f(z) < c . e M<m)M 



x-l 



< c(d+ 1)I 2 I for 





dx™ k 


r x l 


dx\ 


10 


x k 


Jo 


x\-bi--- 6 fc _i 





dx™ 1 


rx\ 


Jb 1 


x l 


Jo 



ldx 



. constant 



for some c € M + (again independent of the choice of z). To conclude that f(z) is of C-type, we 
have to show that we can choose a c G M + such that for all 6 G R we have 



sl - 6i • • • 6 



fc-i 



x * dx\ 



mi 



In 



mi 



"l x bi 



dx 



We show that there is a c G 



such that 



x-l 



< c-el fe l for < x < d and all b G R+: If 



d < 1 we have 



x-l 



< jzg < ■ e' 6 L Finally let d > 1. Since liia^i 



x<"-l 



x-l 



|b| < el fo l there 



are e > and c > such that 



X Cl -l 



x-l 



< ce' b ' for 1 — + It remains to show that 



X<"-1 



x-l 



< cel fe l for x G (0, 1 - e) U (1 +e, d) : 



x"'-l 



x-l 



< - < -e' 6 L As for the proof part concerning 



the exponential growth the integral can be bounded as desired. 



□ 



As a consequence the analytic continuation of S-sums 5' mijm , 2v .. )mfc (&i, 62, ■ . • , 6&; n) with mj G 
N and 6j G K + within the class of C-type functions is uniquely determined by the integral 



mi,r7i2,...,rrifc 



(61,62, • • -,b k ;z). 



If not all bi G M , we can split the integral at whenever the integration interval contains 0. 
For example let us consider the integral representation of Si i 2,2( — 1> 2, — |; n): 



Si, 



2:2 



T,2,--;n 




/o -/II 



24 /"^3 /"^2 



-2 JO 



x 5 (2 + x 4 ) x 3 (1 + x 2 ) (-1 + xi) 



dx\dx2dx^dx^dx 



5- 



Splitting at zero yields 




dx^dx^ 



p0 rx2 




-2 Jx 3 



-l + (-l) n (-xi) n 
x 3 (1 + x 2 ) (-1 + Xl 



1 



+ 



/0 ■/(! 



x 5 (2 + x 4 ) x 3 (1 + x 2 ) 

x 4 />x 3 />x 2 _j _|_ x n 



dx2dx^dxidx5 




dxidx2dx3 



-1 + (-!)' 



-Xl 1 















x 5 (2 + x 4 ) x 3 (1 + x 2 ) (-1 + xi) 



-1 + xi 
dx\dx2dx^,dx4,dx^. 



dxi 
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0.5 




An analytic continuation of (— l) n does not exist (in particular, it is not of C-type), hence we 
cannot continue Si^f - 1> 2, — |; n) analytically. However we can continue the functions 

fe(n) := Si, 2 , 2 f-l,2,-i;2n 

dx^dxs \ f [ f r-t — 1 rdxidx2C?X3 



A) io X b {2 + X 4 ) J J_ 2 Jx 3 J-l X 3 (1 + X 2 ) (-1 + Xi) 

+ ( I I I I — 7^~, r — TT< rdx 2 dx 3 dx 4 dx 5 ] / 1 + ( g rf , r| 

J0 J0 Jo X 5 {2 + X4)X 3 {1 + X 2 ) J + 

+ 11 I I \ — 77i i \ , 1 w , , rdxidx^x^dxidxc, 

/o x 5 (2 + j;4)a;3(l + X2)(-l + xi) 



JO ^0 ^0 



and 

/ (n) := S 1)2 J-l,2,-i;2n + l 



2 




dx^dxz, / / / — — ; — rdxidx 2 dx% 



A) x 5 (2 + x 4 ) / 7-2 7a; 3 7-1 2:3 (1 +x 2 ) (-1 + 



/ / / / — 777- x 77-, rdx 2 dx 3 dx 4 dx 5 / — dxi 

Jo Jo Jo Jo a;5(2 + X4)x 3 (l + X2) y 7 1 -l + a?i 

ft f XS f Xi f X3 f X2 —l+X 2U+l 

+ / / / / / — 77 s 77 — - — 77 — ; -dx\dx 2 dx 3 dx,idx^ 



Jo Jo Jo 



o x 5 (2 + £4)2:3 (1 + x 2 ) (-1 + xi) 



analytically since (—1)™ = 1 for n even and (—1)™ = —1 for n odd. For a graphical illustration 
we refer to Figure [TJ More generally, there is the following result. 

Proposition 6.4. For nii 6 N, tj G M* and z G C consider the integral f(z) = 



Im±,rn 2 ,...,rn k (bi, b 2 , ... , z) from Theorem 2.2 Then f(2z) and f(2z + 1) are of C-type. 



The proof is analogous to Proposition |6.3| Consequently the analytic continuation of even and 
odd 5-sums S mi , m2) ..., mfc (6i, b 2 , . . . , b k ; 2n) and S mi ,m 2 ,...,m k (h, b 2 , . . . , b k ; 2n + 1) with m, G N 
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and bi G M* within the class of C-type functions are uniquely determined with the integral 
representations i mi , m3 ,..., mfc (6i, b 2 , • • • ,b k ;2z) and I mi , m2 ,...,m k (bi,b 2 , ...,£>&; 2* + 1), respectively. 

We remark that in physics applications due to the definiteness of the crossing relations in case 
of the respective process using the light-cone expansion \VS'6\ I14H 1551 E3 EE] n is either even or 
odd and the analytic continuation is carried out from the respective subset. 

7 The Mellin Transform and Inverse Mellin Transform Between 
Multiple Polylogarithms and 5-sums 

For a harmonic polylogarithm h(x), i.e., a generalized polylogarithm with the index set {—1, 0, 1}, 



the Mellin transform (1.1) can be expressed in terms harmonic sums and constants (written 
in terms of harmonic polylogarithms at 1 or equivalently in terms of generalized zeta values). 
Extending the Mellin transform to 

M + (h(x),n) :=M(h{x),n) = [ x n h(x)dx 



M + fM,„) :=M (M,„U 1 I .M dl 
yi + x 1 j yi + x 1 j j i + x 

) Jo l ~ x 

gives a 1-1 correspondence between the harmonic polylogarithms and the harmonic sums via the 
extended Mellin transform and its inverse transform. 

In this Section we generalize this construction and look at the Mellin-transform of generalized 
polylogarithms with indices in K \ (0, 1). It will turn out that these transforms can be expressed 
using a subclass of the S-sums (containing the harmonic sums) together with constants (written 
in terms of generalized polylogarithms evaluated at certain points or equivalently in terms of 
infinite S-sums; see Section [5]). In particular, using the inverse Mellin transform, the 5-sums can 
be expressed in terms of the Mellin transform of generalized polylogarithms. 

In order to establish this connection, we have to deal with functions like f(x) = l/(a — x) 
with a G (0, 1). However, already for this simple case, the Mellin transform is not defined since 
the integral does not converge. Hence we modify the definition of the Mellin transform 

to include these cases, like l/(a 5 ELS follows. 

Let h(x) be a generalized polylogarithm with indices in K\ (0, 1) or h(x) = 1 for x G [0, 1]; let 
a G (0, 00), a\ G (1, 00), 02 G (0, 1]. Then we define the extended and modified Mellin-transform 
as follows: 



M+(h(x),n) : = M(h{x),n) = [ x n 

Jo 



h(x)dx, 

, r , Kx) \ h(x) \ f 1 x n h(x) , 

a + x J \a + x J J a + x 

a\ — x ) \a\ — x ) J a\ — x 

h(x) \ (x n -l)h(a 2 x) ^ f 1 ((^r-l)h(x) 



. = / ^ -'■'^*"' dx= / '—^dx. (7.1) 

\a 2 -x ) J Q l-x J a 2 -x 

In (7.1) both extensions Jq^ 2 \Zx h{a 2 x)dx and ^ x ^^!_ x ^ h(x)dx are equivalent, which can 



be seen using a simple substitution. Since from an algorithmic point of view it is easier to handle 
the second integral, we prefer this representation. From now on we will call the extended and 
modified Mellin transform M + just Mellin transform and we will write M instead of M + . 
For later considerations the following property (see Section [6]) will be helpful. 
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Lemma 7.1. Let m, G R \ (0, 1), a G M* and m e {0, 1} and de/me /(z) = M^ ^^g^ ,, 
T/ien /(z), /(2z) and /(2z + 1) are of C-type. 



7.1 Calculating the Mellin Transform in terms of 5-sums 

Subsequently, we solve the following problem: 

Given m = (mi, . . . , m&) with rrii G R \ (0, 1), a G M* and m G {0, 1}. 

Find an expression F(n) given as a linear combination of S-sums such that for all n G N we have 



Ml ,n 

(a±x) m 



F(n). 



(7.2) 



Note that this construction represents the Mellin transform of weighted generalized polyloga- 
rithms only for positive integer points. However, the uniquely determined analytic continuation 
(among all the C-type functions) for even and odd S-sums describes the Mellin transform in the 



complex plane; further details are given below of (7.5) 



In order to solve the specified problem, we enhance the ideas of [105] for analytic aspects from 
the class of harmonic polylogarithms to the class of generalized polylogarithms. 

Namely, due to the following Lemma we are able to reduce the calculation of the Mellin 
transform of generalized polylogarithms with indices in R\(0, 1) weighted by l/(a+x) or l/(a— x) 
to the calculation of the Mellin transform of generalized polylogarithms with indices in R \ (0, 1) 
which are not weighted, i.e., to the calculation of expressions of the form M(H m (x) ,n) . 



Lemma 7.2. Forn G N, m = (mi, m) 
a\ G (l,oo), a2 G (0,1), we have 



(mi, m2, • • • , rrik) with mj G K\(0, 1) and a G (0, oo) 



M 
M 



1 



a + x 
1 

a — x 



, n 



, n 



a + x 

N[['5=M,n 
1 — x 

• | 

a\ — x 



(-a)»(Si(-i;n)+H_ (l)), 

(_ a )»( Sl (-J;„) 



if < a < 1 



Si(-i;oo)), ifl>a 

if < a < 1 
Si(^)+Si(i;oo)), ifl<« 



-n M(H_ a , m (x) , n - 1) + H_ a , m (l) 

-n M(Hi jm (x),n-l), 

-n M(H ai , m (x) , n - 1) + H ai , m (l) . 



a 2 - x 



n 



-H mij m(l) S] 



-:n 



En 
i=l 



' 1 Vlllf H " W i 



«2 

'2 / \ l-^ 



-H 0)7fi (l)Si(i;n)+Er=i^ 

(— 

v^n V_a 



M 



ggW 



-H mi ^(l)Si^;nJ+Er=i 
where the arising constants on the right hand side are finite. 
Proof. For a > we get 

i \ r 1 x n , r 1 x n - {-a) n 



i Vm 



mi —x 1 



if mi < 

if mi = 

if mi = 1, 
if mi > 1, 



M 



a + x 



, n 



o a + x 
o 



-dx 



a + x 



dx+ I —dx 

a + x 



n-l 



i=0 ^ ' 
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n— 1 -. 

(-«)"- 1 E(Z iP (i T T) + (-°}" H -«M 

(-a)"Si +(-a)"H_ a (l). 



If a > 1 this is equal to (— a) n Si(-^; n) - Si(— \; oo) 
For < a < 1 we get 



"(—•») = -J* x — = / 1 «- 1 E( 1 )'* 

\a-x y a n J a-x a n J 

n—l , ,i , n—l 



- V - / x*dx = - V 

a t J al Jo a^a^l + i) 

-Si(-;n 



For a > 1 we get 



l \ f 1 x™ x n - a n Z" 1 1 

M( ,n) = / dx = / dx + a n dx 

,a-x J J Q a-x J a-x J a-x 

-a n Si ( -; n ) + a n H a (l) = -a n Si ( - ; n j + a n Si ( -; oo 



For a > we get 



n + 1 n + 1 \ a + x / 



M _ H4x)^ + . _( n + i) M (H_ aim (x),n) + H_ a , m (l). 



Hence we get 



n a + x 

Similarly we get 

"1 i / /•£ „n+l 



[ x n H 1>m (x) dx = — l — lim f e" +1 H ljm (e) - f ^ -H m (x) dx + H 1 , m (e)') 

7o n + i v y 1 - X 'J 



1 /•! _ 1 

— — r / —. R m (x)dx 

n + 1 J 1- x 

1 U(^,n + 1 ). 



n + 1 V 1 — x 

And hence 



( JT^U ■ 1 : - -(, H UirCfl,,,,,,!,-).,;). 



For ai € (1, oo) we get 

H ai , m (l) 1 f l x n+1 



f x"H ai , m (x)dx = H ^ 1) -^ T / 
Jo n + 1 n + 1 J (ai - x) 



H m (x) dx 



23 



Hence we get 



M 



_ 3a ! =W_^ M f5=M 1 , + i y 

n + 1 n + 1 \ai — x / 

f H ^,n+l) = -(n+l)M(H ai , m (x),n)+H ai , m (l). 
y 01 x y 



For a2 £ (0, 1) and mi < we get 



M( Hmi ' m(x) ,n 

(12 — X 



1 f 1 (x n - a£)H mi , m (x) 



4/ 

J 2 ./0 



<ix 



a2 — x 



— H mi m (x) 



Z" 1 H m (x) 

7 |mi|+x^ i 



-dx 



Tl M\Q / 1 ^ J V a2 / / X*H m X 

-H mi , m 1 Si — ;n + > — : — / | j- — dx 

\a 2 J f- { % Jo mi +35 



1 



-H m i,m(l)Sif — ;nj 

i=i 



<w M ( H m (x) 



|mi| + x : 



and 



a2 - x 



a 2 - x 



1 Z" 1 (x n 



(x n - a2)H , m (x) 



dx 



a2 — x 



-H , m (x) ^ 



(J 2 



1=1 



= -H 0>m (l) Si( — ;r; 



+ 



"2 



1 H m (x) \ -> 

<> x h % 



-dx 



i=i 



a 2 J f 1 X J H m (x) 



dx 



* Jo x 



_ _ Ho , m (l) Sl |I;„) + ^4I M (5=W, !l , 



1 Um /'(x'-^..^^ 



a, e-H 



a 2 - x 



= J™ I -Hi, m (x) ^ 



lim [ -Hi jm (e) 



r H m (x) ^ 
Jo 

ji.iyti r ^H m (x) 



dx 



= lim^-Hi, ra (e)Si^— ;nj +Si 
+ y- KL r (g ~ l)H m (x) 



Hi, m (e) 



dx 



24 



= £ 

For 02 £ (0, 1) and mi > we get 



-M 



H m (x) 
1-x 



M 



H 



mi,mV 



a 2 - x ' / a£ Jo a 2 - x 



-H mijm (x) 



= -H 



— — H mi>m (l)Si 



Z" 1 H m (x) (a 2 ) 
Jo m l ~ x ~[ i 

V a 2 / « Jo m l~X 

n 



dx 



(fx 



\mi — x ' 



□ 



To calculate a non- weighted Mellin transform M(H m (x) , n) we proceed by recursion. First let 
us state the base cases, i.e., the Mellin transforms of generalized polylogarithms with depth 1. 

Lemma 7.3. For n G N, a G R \ (0, 1) we have 

{ (^w(l + « n+ V + l)Si(i;n) 

-(a" +1 -l)(n + l)Si(i;oo)), ifa<-l 



M(H a (x),n) = { 



^(l + a^n + lJSiQjn) 

+ (a" +1 -l)(n + l)H a (l)), 



(n+l) 2 ' 

^(l + (n + l)Si(n)), 



if - 1< a < 
if a = 
if a = 1 



(^(i + ^n + lJSiQjn) 

_( a n+l_l)( n + l) Sl (I;oo)), ifa>l 

where the arising constants are finite. 
Proof. First let a < —1. By integration by parts we get: 

M(H a (x),n) = / x n R a (x) dx = — -H a ( 
Jo n + l 

- h.<L> -J-(f ^'-^'rix + „•» /' ' ^ 

1 n + l |a|+x Jo |a|+x / 



o Jo 



1 x n+l 



1 



n + 1 lal + x 



-<ix 



n + 
Ho(l) 



(1 _ a n+l) + 

n + 1 v y n + 



l « 



Ha(l) 



(1-a 



n+1^ 



n + l 

-SiQ;oc 



n+l ^ aHi + l) 

t=0 



n + l 



-(1 -a n+l ) 



,71+1 



n + l 



Si -;n + l 
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n + l V n + 1 



+ a n+1 S 1 [-;n)-(a n+1 -l)S 1 [-; 



; oo 



For — 1 < a < we obtain: 
M(H a (x),n) = 



1 x n+l 

x n H a (x) dx = — — 
o n + l 



h„(x; 



1 f 1 x n+1 1 
J n+l\a\- 



1 



Hg(l) 

n + 1 n + 1 

Ha(l) 



1 ^.n+l _ a n +l 



dx + a n+1 
a| + x ,/o l a l + x 



+ x 
1 



dx 



a" +1 )-— I 1 a n Y- 
n + 1 n + 1 7 ^ a* 

Ha(1) (l a"+M ^ V 1 



■(1 



n + 

Ha(l) 

n + 1 
-1 



(l_ a n+i)_ ^— -Sil -;n + l 



i=0 
n+l / j 



n + l V a 



(n+l) 5 



l + a n+i (n + l)Si( -;n) + (a n+1 - l)(n + l)H a (l) ) . 



For a = it follows that 

M(H (x) , n) 



rl x n+l 1 f 1 X n+1 1 

/ x n H (x) dx = H (x) - / dx 

Jo n + 1 Jo n + l x 

f 1 x n _ x n+1 
~J nTT ~ ~{n~+Tf 



(n + l) 2 ' 



For a = 1 we get: 



f x n Hi(x)dx = lim f x n H!(x)dx 
Jo t-* 1 Jo 



„n+l 



lim 

e-*-i V n + 



T Hi(x) 



„n+l 



n 



n + 



o Jo (n + l)(l-x) 

1 / f £ r n+1 — 1 \ 
— lim e n+1 Hi(e)- / xdx - Hi(e) J 

^-j- ^lim(e n+1 - l)Hi(e) + lim jf x<da ^ 



i i i 

XT E -ri = r^TY2 (! + (» + !) s i (")) • 

n + l^-^i + i (n + lr 

i=0 

Finally for a > 1 we conclude that 

M(H (x),n) = / x"H a (x)dx = — -H a (x) -/ — — 

Jo n + 1 Jo n + la - x 



dx 



Hq(l) 1 

n+l n+l 

n + 



/ r 1 x n + l a^ dx + a „ +1 r 1 _J_ d \ 

\Jo a - x Jo a - x J 



a — x 

i n Ji 



1 n + l Jo ^a* 
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Ha(l) 



n + 



i(l-a« +1 ) + ^-y — - 



(i + 1) 

^(l-a-) + ^ Sl (l;n + l 
n + 1 n + 1 V a 



— f— + a" +1 S 1 f 1 ;nV( 
n+l \n+l \a J V 



a n+1 - l)Si( ^;oo 



The higher depth results for M(H m (x) , n) can now be obtained by recursion: 
Lemma 7.4. For n G N, a G K \ (0, 1) and m € (R \ (0, l)) fe , 



M(H a , m (x),n) 



f (l~a" +1 )H a ^l) _ V n MH^),«) : f _ ^ n 

n+l n+l l^i=0 a i ' 11 « ^ u 

^-^(H^n), if O = 

^5X M(H m (x),n), ifo = l 

(l-a^H^l) _a-_ v „ MH^),t) ., , 
n+l n+l Z^j=0 a* ' 11 " ^ 1 

where the arising constants are finite. 

Proof. We get the following results by integration by parts. For a < we get: 

M(H a , m (x) ,n) = x n H flim (x) dx = — — H ajm (x) 
Jo n+l 

H a , m (l) 1 



1 ^ x n+1 H m ,fx) 1 



n+l n+l 



1 x n+1 - a n+1 



j n + l \a\ + x 
1 H m (x) 



-dx 



H m (x)dx + a n+1 / 

JO 



H a , m (l) _ a n+l\ _ 

n+l 1 ; n + 



|a| + x 

"' (a) 

a* 



a + x 



dx 



_ H a;m (l) ^ a n + ls ) _ Q " 
n+l 1 ; n + 



i=0 



For a = we get: 



f 1 X n R ,r 

Jo 



,{x) dx = 



x n+1 

^TT Ho ' m(x: 

H ,m(l) I 



1 rl 



n+l n+l 



o Jo n + 1 
M(H m (x) , n) . 



H m (x) dx 



For a = 1 we it follows that 



/ x n Hi m (x) dx = lim / x ra Hi m (x) dx 
Jo ^ Jo 



-H m (x) dx 



/ x n+l e r 

e^i\n + l 7 (n + l)(l-x) 

1 / t' e r n+l — 1 \ 

= — — lim e" +1 H lim (e) - / — H m (x) dx - H ljfn (e) 

n + 1 e->i V J 1-x J 

= —J— f lim(e" +1 - l)Hi, m (e) + lim f V x l H m (x)dx ] 

n + i^i ^i7 ^ y 



x l H m (x)dx 



□ 



27 



1 n 

i=0 



For a > 1 we conclude that 



„n+l 



M(H am (x),n) = / i ra H am (g) rfx = ~ H a m (x) 
Jo n + 1 



1 /•! „n+l 



JO 



s n+1 H m (s) 1 



n + 1 a — x 



-dx 



H a , m (l) 1 



n + 1 n + 1 



1 x n+l - a n+l 



a — x 



H m (x) dx + a r ' 



1 f'Umix) 



dx 



a — x 



H a , m (i) (1 _ „+i } ^_ /vv 

J{J i=0 

^(1 - « n+1 ) + E ^M( H -^) , • 

1 n + 1 a* 

i=0 



1 n 



^H m (x) 



n + 

H tt)TO ( l) 
n + 



□ 



Summarizing, using Lemma 7.2 together with Lemma 7.3 and Lemma |7.4| we are able to 
calculate the Mellin transform of generalized poly logarithms with indices in R\ (0, 1). In addition, 
the polylogarithms can be weighted by l/(a + x) or l/(o — x) for a£l. These Mellin transforms 
can be expressed using S*-sums. 

As an example we compute the Mellin transform of 

x + 2 



Using Lemma 7.2 we can write 

'H _ 2 (x 



M 



x + 2 



n 



-n • M(H_ 2A _ 2 (x) , n - 1) + H_ 2)0 ,_ 2 (1) 



Now we apply Lemma 7.4 which leads to 

1 - 



M(H_ 2 ,o,-2(x),n-l) 



n 



2)n H m ( - 2)n_1 

H_2,0,-2(.J-j — 



n 



n-1 

£ 

i=0 



M(Ho,- 2 (a;),i) 



(-2)* 



Again using Lemma 7.4 we find 

M(H ,_ 2 (x) , i) = J-Ho.-aCl) -4tE M ( H -2^) , • 
i + l i + l ~ i 

Finally applying Lemma |7.3| we get 

M(H_ 2 (x),i) = 



i=0 



1 + (-2)*_(i + l)Si(-^; Q - ((-2)* 1 - l)(i + l)Si(-^; oo)) 

(i + l) 2 



Putting this together yields 
'H ,- 2 (x) 



M 



x + 2 



,n\ =(-2)"H_ 2 n,- 2 (l) 



+ (-2r- 1 E^(^ 2 (D + 

i=0 



l+(-2)W(l+QS 1 (-|;»)-(-l+(-2) 1 +')(l+i)Si(-l;cx,) 
(1+i) 2 



)■ 
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To this end, transforming this expression to S-sums we arrive at 
'Ho,- 2 (x) 



M 



x + 2 



, n 



(-l)»(2) n H_ 2 ,o,- 2 (l) + Ho,_ 2 (l)Si (-i; n) - S 2 (n)Si (-§; oo) 



+ Si (-5500) S 2 (-|;n) + S 2 ,i (l, -fjn) 



(7.4) 



Looking at the presented algorithm the representation F(n) in (7.2) will be of the form 
F(n) = (-l) n La5/o(n) + • • • + r d a n d f d {n)] + L+K^+i/d+lfa) + ■ • • + r e a n J e {n)\ (7.5) 



with d < e, r, G K, aj G M + and where each /j(n) stands for a S-sum. 



Consider the C-type function f(z) := M(t^ 



H m {x) 



see Lemma 







We remark that the 



right hand side of (7.2) for even and odd n has also an analytic continuation by replacing the S- 
sums with the integral representation of Theorem 2.2 In particular the function is of C-type due 
to Proposition |6.4| and Lemma 6.2 Because of Carlson's Theorem it follows that the constructed 
function agrees with f(z). In other words, the representation of the Mellin transform with S-sums 
for integer values gives the representation in the complex plane. 

7.2 Calculating the Inverse Mellin Transform in terms of generalized polylog- 
arithms 

Subsequently, we consider a subset of the generalized harmonic polylogarithms, namely the set 



H 



H 



mi,m2,...,m,f 1 



(X) 





1 






mi G R, 


< 1 and 


m h 


Win. , 
Jl — 1 



< 1 for 2 < i < k where 



{rrij! , rrij 2 , . . . , rrij k } is the set of nonzero indices, and j u < j v if u < v 



and call the elements ^-polylogarithms; note that m; G R\(— 1, 1) with mi 7^ 0. For this subclass, 
it turns out that the Mellin transform of properly weighted generalized polylogarithms can be 
expressed using the S'-sums from the set 



S := < S ai) a 2 ,...,a h (bi,b2, ...,b k ;n) 



(H G IT for 1 < t < k; b x G W 
bi G [-1,1] \ {0} for 2 < i < k }; 



(7.6) 



the elements of this set are also called S-sums. 

Restricting to this class, we solve the following problem: 
Given an S-sum S m (fa; n). 

Find a linear combination F(n) of Mellin transforms of weighted generalized polylogarithms from 
H such that for all n G N we have 

S m (6;n) = F{n). (7.7) 

In other words, there is a certain bijection between the set of properly weighted H— polylogarithms 
and S-sums via the Mellin-transform and its inverse Mellin transform. 

In order to present the algorithm for the inverse Mellin transform, we define the following 
order on S'-sums. Let S mi (&i;n) and S m2 (b2',n) be S-sums with weights w\, W2 and depths d\ 
and d 2 , respectively. Then 

S mi (&i; n) -< S m2 (b 2 ; n) , if w\ < w 2 , or (w\ = w 2 and d\ < d 2 ). 

We say that a S-sum s\ is more complicated than a S-sum s 2 if s 2 -< s%. For a set of S-sums we 
call a S-sum most complicated if it is a largest element with respect to -<. 

Then the inverse Mellin transform from S to H relies on the following main properties: 
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Lemma 7.5. 

1. In the Mellin transform of a generalized polylogarithm H m (a?) G H with m = 
(roi,m2, • • • , rrifc) weighted by l/(c — ») or l/(c + x) (c G M*j i/iere is or% one most com- 
plicated S-sum S a (b; n) G S with a= (ai, . . . , aj) and b = (&i, . . . , 6j) ,ie., 

M /Hm^) A =p . Sa(b;n) + t (7 . 8) 

\ C ± X / 

where p£l* and i consists of a linear combination of S-sums over R i/iai are smaller than 
S a (b; n). Note that ifH m (x) is of weight w, then S a (fa; n) is of weight w + 1. In particular, 
if mk 7^ 1 fien a; 7^ 1 or 6; / 1; i/zese properties will be used later in Section\Tl\ 

2. Given S a (b; n) G S with a = (ai, . . . , a;) and 6 = (61, ... , 6;), i/iere is an algorithm (see 
Algorithms 3]) that computes an H m (x) G H with m= (mi, ra2, • . . , mi-), d£ {—1,1}, and 

c G M* such that S a (b;n) is the most complicated S-sum that occurs in M^ ^^l , 71J . Note 
that if S a (6;n) is 0/ weight uu, then H m (x) is of weight w — 1. 7n particular, if a\ 7^ 1 or 
bi ^= 1 then 7^ 1. Moreover, if bi G [—1, 1] loit/t 6j 7^ /or a// i f/ien c > 1; again these 
properties will be used later in Section [71} 

Note that parts of these properties have been exploited already in |105| to construct the 
inverse Mellin transform between harmonic sums and harmonic poly logarithms. The rigorous 
proofs for this special case given in [1] carry over to this more general case and are omitted here. 

The computation of the inverse Mellin transform of a linear combination of S-sums now is 
straightforward (compare [H 105J): 

• Locate the most complicated S-sum. 

• Exploit property (2) from above and compute an ^-polylogarithm weighted by l/(c — x) 
or l/(c + x) such that the S-sum is the most complicated S-sum in the Mellin transform of 
this weighted generalized polylogarithm. 

• Add it and subtract it with an appropriate factor (see the next step). 

• Perform the Mellin transform to the subtracted version. This will cancel the original S-sum 
(by choosing the appropriate factor). 

• Repeat the above steps until there are no more S-sums. 

• Let c be the remaining constant term; replace it by M _1 (c), or equivalently, multiply c by 
6(1 — x) where S(x) stands for the Dirac-<5-distribution [138J. 

In the following we will illustrate the computation of the inverse Mellin transformation by 
means of an example. We consider the sum S 2j i (l,— |;n). It is the most complicated S-sum in 
the Mellin transformation of H ° x ^_ 2 2 ^ which can be computed by Algorithm 3 in [6]. So in a first 
step we get 

S 2>1 (l,-i;n) = S 2il (1,-I;n) - {-\) n m(^%#, n) + (4)» M (5*#,n) . 
Calculating the first of the two Mellin transforms (see (7.4)) we get 

S 2jl (l,-|;n) = H) n M(^^,n) +S 2 (n)Si (-i;oo) - S 2 (-i;n) Si (-|;oo) 
+ Si (-\;n) S 2 (-|;oo) + S3 (-5500) - Si, 2 (-5, 1; 00) . 
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In the right hand side of this expression we find a linear combination of the three 5-sums 
S2 ( — i 82(^1) an d Si (— 2> n )) which we want to express using the Mellin transformation of 
generalized poly logarithms. The most complicated 5-sums from these three is S2 (— \ ; nj . In par- 
ticular, it is the most complicated S-sum in the Mellin transformation of ■ Hence following 
the above method we rewrite S 2i i (l 

S 2 



+ 



^;n) to the expression 

H;»)-(-§)"m(^)h, „ ... Vx+2 , 

■i; 00) + Si (-|;n) S 2 



1 >) n M(' I ^h#,n)+S2(n)S 



2' 



00 



x + 2 
- Si 2 (-1,1; 00) . 



Si (-5500) 
(-5; 00) 



Performing the first of these occurring Mellin transformations yields 

'H (X) 



i)"M(' H °'- 2(x) 



, n 



+ 



2) 



2 , „^ x + 2 , v 2 , SiH;oc)M Vx + 2 ,, 7 2 , 

-Si,2 (-5, l;oo) + Si (-|;n) S 2 (-|;oo) - Si (— §;oo) S 2 (-|;oo) + S 3 (-5; 00) 
Continuing this process for the remaining 5-sums we arrive at 



+ S 2 (n)Si (-|;oo) 



S 2 ,i (l,-i;n)= Sl (-J;oo)M(^ 

'H (x 



+ (-1) 



n f 1\« 
2/ 



Si 



;oo M 



x + 2 



s 3 (- 
+ s 2 



2 - 



;oo) - Si, 2 (-2,1; 00) 



; ;oo) M 



1 



x + 2' 



, n 



+ M 



H 



0,-2 



(x) 



x + 2 



(7.9) 



We remark that there is a second way to compute the inverse Mellin transform of an 5-sum by 
using the integral representation of Theorem |2.2| Repeated suitable integration by parts of the 
given integral representation provides the inverse Mellin transform in terms of integrals that can 
be expressed in terms of generalized polylogarithms. 



Looking at the presented algorithm the representation F(n) in (7.7) will be of the form 



F(n) = (-l) n [r a%f (n) + ■■■ + r d a n d f d (n)] + [r d+1 a n d+1 f d+1 (n) + ■■■ + r e a n J e (n)] (7.10) 



with d < e, Ti G 



and where each fi(n) stands for a Mellin transform of weighted 



harmonic polylogarithm from H. 



We remark that the right hand side of (7.7) for even and odd n gives rise to an analytic 



continuation. In particular the function is of C-type due to Lemma 7.1 and Lemma 6.2 Because 



of the generalized Carlson Theorem it follows that the obtained function agrees with the integral 
representation of the even, resp., odd 5-sum given in Theorem |2.2| 

In many QCD applications one starts with a complicated function f(x) in x-space where 
one knows that it can be expressed in terms of a linear combination of weighted generalized 
polylogarithms (or related classes of them) over real numbers. It is possible to choose either 
n as an even or odd positive integer, depending on the physics problem. In some cases the 
corresponding other choice may even yield a vanishing result due to a symmetry inherent in the 
process. Then one tactic is to perform the simplification on M(/(x),n) (e.g., with the symbolic 



summation packages [U [221 EB E] ) an d to express it in terms of S'-sums. Note that by Lemma 7.1 
M(/(x),n) is of C-type. Finally, applying the inverse Mellin transform to each of these 5-sums 



and combining these calculations, one obtains an expression F{n) in the format (7.7). Usually, the 
alternating sign and the powers a" vanish. Since also the output function F{n) can be analytically 
continued, the continuation is of C-type and M(/(x),n) = F(n) for almost all integer n, their 
analytic continuations agree by Carlson's Theorem. Now drop the Mellin operator from F(n) 
which gives a linear combination h{x) of generalized polylogarithms such that M(h(x),n) = F(n). 
This finally yields that f(x) = h(x). 
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8 Differentiation of 5-Sums 



As worked out in Section [6] the even and odd S-sums can be uniquely analytically continued 
within the class of C-type functions. This allows one to consider differentiation with respect to n. 
As it turns out the obtained expressions are again integral representations of S-sum expressions. 
This will make it possible to define the differentiation of S-sums in an analytically consistent way. 

We present two strategies to compute the differentiation of S'-sums. The first approach follows 
the ideas from \27\ [26] for harmonic sums: compute the Mellin transform, differentiate the arising 
integrals (which can be analytically continued in n), and apply afterwards the inverse Mellin 
transform to bring these expressions back to S'-sums. This method is restricted to S-sums (for 



a definition see (7.6)) using our current technologies. Finally, we will present a second approach 
using directly the integral representation of S-sums that will work in general. 

In both approaches one can start with an even S-sum S a (b; 2n) or odd S-sum S a (b; 2n + 1) 
(where the uniqueness property of the analytic continuation is given). For convenience, one can 
also work with the usual S-sum and obtains an output of the form 

S a (b;n) = (-l) n F(n) + G(n), 

which encodes both cases: Setting n — > 2n, resp. n — > 2n + 1, yields the differentiation (uniquely 
determined among the C-type functions) for the even resp. odd case. 

8.1 Differentiation by means of the Mellin transform and its inverse 



Using the inverse Mellin transform we express an S-sum to the form (7.10) where the fi(n) stands 
for an expression of the form M ^ , nj with H m (x) 6 H. Note that from now on the possibly 
occurring ( — l) n in (7.10) remains untouched and we deal only with the remaining objects which 
can be analytically continued. Next we differentiate these summands. For "+" and k £ IR + this 
leads to: 

yw.) f v »" i -'\ 

an \k + x J an J k + x J k + x 

X {Ho,mi,...,mi{ x ) ~\~ Hmi,0,m2...,mi (^) + " ' ' + H miji .^m^O^)) 

o k + x 



M( H , T n 1 ,...,m,(a:) + . . . + M ( Hmi,.,mM „ 

k + x J \ k + x 



For "— " and k 6 (1, oo) this leads to: 



— M , n =— / — dx= dx 

an \ k — x J an J k — x J k — x 

^ Ho,mi,...,m; (2") H ) 7j, 1) o,rri<2... ) m; ~\~ ' ' ' ~\~ Hmi,...,mj,o(^) 

o k — x 

= M / H , m m ,(x) \ + + M / H m mi , (x) n 
\ k — x ) \ k — x 

Finally, for "-" and k £ (0, 1] this leads to: 

d M ( Km(x) \ = d r 1 (£-l)R m (x) dx= f 1 g(Ho(x)-Ho(fc))H m (x) ^ 
dn \ k — x ' / dn J Q k — x J Q k — x 

1 (g-w^j. _ H 0( *) /' (g ; 1)H - (I) ^ 

o k — x Jo k — x 
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f 1 H (f ) H m (s) , . 

I / V ' dx. 8.1 

In the following lemmas we will see that ~^~^~—^-dx is finite and in other words, the last 
integral in (8.1 ) is well defined. Moreover, we will see how we can evaluate this constant in terms 
of generalized polylogarithms. Hence in all three cases we eventually arrive at Mellin transforms 
of generalized polylogarithms. To this end, we perform the Mellin transform of these expressions 
and arrive at a final expression in terms of 5-sums together with infinite S-sums. 

Lemma 8.1. For a generalized poll/logarithm H m (x) with m = (mi, m-2, . . . , mi), rrii £ M\ (0, 1] 
and i 6 N the integral 

nl H ( x)* Km(x) dx 
1 — X 

exists and can be expressed in terms of generalized polylogarithms. 
Proof. Using integration by parts we get 



1 — X a-S>0+ J n 1 — X 



lim lim [ ^(x) 4 Hi m (x) & - i T ^'"^ rix 

6->l- a^0+ \ a J a X 

I' 

JO 



. Z 1 H (x) t_1 H lim (x) , 
/ ax. 

x 



After expanding the product Ho(x)* 1 Hi )TO (x) and applying the integral we end up in generalized 
polylogarithms at one with leading zero. Hence the integral is finite. □ 

Lemma 8.2. For k G (0, 1) and a = (a%, <Z2, ■ ■ ■ , ai), m;£K \ (0, 1] we have that 

- 1 Ho(f)H a (x) ^ 
— x 

is finite and can be expressed in terms of generalized polylogarithms. 

Proof. If H a (x) has trailing zeroes, i.e., ai = 0, we first extract them and get a linear combination 
of expressions of the form 



1 H (f)H (xrH m (x) 
k — x 

with w GN. Now let us look at such an expression: 



dx 



1 H (f)H (xrH m (x)^ ri K (x)Ko(kx) w H m (kx) kdx 



k x _/n k kx 

| H (x) (H (x) + HotAO^Hmj. m (x) 

k ' ' k dx 

1 — x 

| H^ ^(x) Eto (T)Ho(^) i+1 Ko(k) w - 1 

k v "' k 

1 — X 

' W \ _, |i% rn L ( X )R (xy + ' 1 

i=0 



r/.r 
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riHmj. ^(x)H (x) l+1 /•iHm 1 ^(x)H (x) ,+1 

A = — fc — fc dx+ / — fc — k dx 

Jo 1-x J i l-x 

•iH™, mt(s)Ho(x) i+1 fE-iHmi m l (x + l)H (x + l) i+1 
— fc fc dx— / — fc k dx . 



o x Jo x 



B C 

The integral -B is finite and expressible in terms of generalized polylogarithms due to the prev ious 



3.1 



and 



lemma. After expanding the product Hmi ™±(x + 1) Ho(x + l) 4 , applying Lemma 

k fc 

applying the integral we end up in generalized polylogarithms at one and at -r — 1 which are all 
finite. Hence C is finite. Summarizing 



L 



io k — x 

is finite, since we can write it as a sum of finite integrals. □ 
We demonstrate this strategy by differentiating S2(2;n). After calculating 



2 % J Jo 2 X 



we differentiate with respect to n and obtain 

d f 1 (2 n x n - l)Ho(x) , f 1 Tx n (H (x) + H (2)) H (x) 



dx = 5 dx 



|-x Jo |-x 

+ /' H °' 2 ^ H ° (I) & 

Jo 2~ X 

2m (^,„Uh„(2)m('^,„) + ^Ho(2,)Ho(x )[fa 



2 / V 2 / " 2 *^ 



= -2S 3 (2; n) + H (2) S 2 (2; n) + H ,o,-i(l) + 2H 0j o,i(l) + H ,i,-i(l) 
8.2 Differentiation based on the nested integral representation 



Using the integral representation from Theorem 2.2 the differentiation can be extended from S- 
sums to S-sums. We demonstrate the underlying ideas by the differentiation of Si^iQ, 2, |; n). 
We start with the iterated integral representation 

i 

f3 f x [ v 1 f z w n — 1 

S w (i,2,i;n) =J^ ^ ^ \ ^ — j J 1 ^-dwdzdydx 



2 2 

si r 1 i r y i r- ,(•■■ 



x - 1 J x y Jo l-zjl w-l 



dwdzdydx 



given by Theorem 2.2 Differentiation with respect to n yields 

rl i r y i r w n n (w) 



o x - 1 J x y J 1 - z J\ w-l 



dwdzdydx =: A. 



2 
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We want to rewrite A in terms of S-sums at n and finite S'-sums at oo (or finite generalized 
polylogarithms) . Therefore we first rewrite A in the following form: 

rk i r 1 i rv i f z w n -l 
A = J Q - y j Q —M w)dwdzdydx + 

S v ' 

S: = 

/ — *— [ -- [ y~ — / H °( w ) dwdzdydx . 
Jo x-1 J x V Jo z Jl w-1 

v v ' 

C: = 

Let us first consider 5: we split the integral at the zeroes of the denominators and get 




— Bi + £?2 + B3. 

Starting from the inner integral and integrating integral by integral leads to 




+S 2 (|;n) - S 2 (z;n) \ dzdydx 



= -H 3 (l)H ,i,o(l)S 1 (i;n) +H (2)H 3j0j 3(l)Si (±; n) - H (3)H 3 3 (l)Si (§;n) 

- H o,3,o,| (1)Sl _ 2H 3,o A | (1)Sl ^ ;n ) _ H o(2)H 3 (l)S M (i, l;n) 

+Ho(3)H3(l)Si, 2 (|,2;n)+Ho,3(l)Si,2(|,2;n) -2H 3 (l)Si, 3 (|,l;n) 
+2H 3 (l)S li3 (|, 2; n) - H 3 (l)S 2 , 2 (i, 1; n) + H 3 (1)S 2 , 2 (|, 2; n) 
-Ho(3)Si, 2 ,i(i,2,|;n) -Si, 2 , 2 (|,2,|;n) -2Si, 3 ,i(|,2,|;n) 
— S 2 , 2,i( 2 , 2, |; n) . 

Applying the same strategy to B2 and B% leads to 

B = -H 3 (l)Ho,i, (2)Si (|;n) +^(2^3(1)8! (§; n) - H (3)H 3 3 (l)Si (§;n) 

- H o,3,o,| (1)Sl ^ ;n ) _ 2H 3,oa1 (1)Si (|;«) +Ho(3)H 3 (l)Si, 2 (|,2;n) 
+H ,3(l)Si, 2 (|, 2; n) - H (3)Si, 2 ,i(|, 2, §; n) - S lj2 , 2 (i, 2, §; n) 
-2Si,3,i(|,2,|;n)-S 2l 2,i(|,2,|;n). 

Let us now consider C: again we split the integral at the zeroes of the denominators and get: 

c = r j- r - 1 - r '-^ r -^i dwdzdydx+ 

Jo x - 1 Jx V Jo z- ± J z w-1 

n 1 r 1 1 n 1 p H («;) J , , , 
Jo —iJ^yJo —J z -—^^^ 

lo^!\- l y!lz^\ II ^*°***>*° 
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— C\ + C*2 + C3. 

Starting from the inner integral and integrating integral by integral leads to 



^ — J x y Jo z 2 

1 rh -Hi(y)H li0 (|)+Hx )li0 (y) 
— / ? = dydx 

^ — 1 Jx y 

Hi.o(§) (Hid) - H - H o|>10 (i) + H o i 10 W 

dx 

x-l 

H (2)H 2 (l)H 3 (l)H ,i(l) + H 3 (l)H 0)2 (l)H ,i(l) - H (2)H 3 (l)Ho,i, 2 (l) 
-Hq(2)H 2 (1)H n 3 (1) - H 0)2 (1)H n 3 (1) - 2H 3 (l)H 0> o,i l2 (l) - H 3 (l)H ,i, ,2(l) 



+H (3)H 3 (1) + H 3 (1) + 2H 3 (1) + H 3 (1). 

n ' 3,0,§,3 V ; 0,3,0,§,3 V ' 3,0,0,§,3 V ; 3,0,|,0,3 V ' 

Applying the same strategy to C2 and C3 leads to 

C = H 3 (l)Ho(2) 2 H ,-i(l) + H 2 (l)H 3 (l)Ho(2)H ,i(l) + H 3 (l)H (2)H_ 1>0 ,i(l) 
+H 3 (l)H (2)H ,_i, 1 (l) - H 3 (l)Ho(2)Ho,i,2(l) - H 2 (1)H (2)H 3 (1) 

3,U, 2 

+H 3 (l)Ho,i(l)Ho, 2 (l) - H , 2 (1)H 3(1) - H 3 (l)H_ lj0 ,i,_i(l) - 2H 3 (1)H 0A1 , 2 (1) 



0,3,0,^,3 V ' 3,0,0, 5, 3 V ' 3,0, |, 0,3' 



-H 3 (l)Ho il)0>2 (l) + H (3)H Q 3 3 (1) +H n , n a ,(1) + 2H, n „ a ,(1) + H, n 3 n ,(1). 

Adding B and C leads to the final result: 

aSi, 2 ,i(f,2,§;n) _ 



dn 

Ho(2)H 3j0j s(1)Si (|; n) - H 3 (l)H ,i 1 d(2)S 1 (§; n) - Hq(3)H 3 Q 3 (l)Si (|;n) 

- \ itQt dX)*i (|;») - 2H 3 0) 3(1)S! + Ho(3)Ha(l)Si, 2 (|, 2; n) 
+ H , 3 (l)S 1 , 2 (i,2;n) -Ho(3)Si, 2i i(|,2, |;n) _ Si, 2i2 (|,2, |;n) 

- 2Si, 3 ,i(|, 2, |; n) - S 2)2ll (i, 2, §; n) + H 3 (l)Ho(2) 2 H ,_i(l) 

+ H 2 (l)H 3 (l)H (2)H ,i(l) + H 3 (l)H (2)H_ li0 ,i(l) + H 3 (l)H (2)H ,-i,i(l) 

- H3(l)Ho(2)Ho,i,d(l) - H 2 (1)H (2)H 3 0)3 (1) + H 3 (l)H oa (l)H , 2 (l) - H 0>2 (1)H 3 >Q 3 (1) 

- H3(l)H_ li0 ,i,_j(l) - 2H 3 (l)H ,o,i,2(l) - H3(l)Ho,i,o^l) 



(8.2) 



+ H (3)H 3 (1) + H 3 (1)+2H 3 (1)+H 3 (1). 

y ' 3,0,§,3 V ; 0,3,0, §,3 V ; 3,0,0, §,3 V ; 3,0,f,0,3 v ' 

As worked out in [6] the proposed strategy works in general for S-sums. 

9 Relations between S'-Sums 

iS-sums obey stuflie relations due to their quasi-shuffle algebra, differentiation relations, and 
generalized argument relations. Subsequently, the available ideas for harmonic sums |21|, [2?1 [26] 
and cyclotomic harmonic sums [5] are extended and explored for S'-sums. 

9.1 Quasi-Shuffle or Stuffle Relations 

Based on [7T] it has been worked out in [90J that Z-sums 

Z ai ,...,a k (xi,...,x k ;n)= V _i....-A. (9.1) 

1 * If 



36 



with their quasi-shuffle product form a quasi-shuffle algebra. Since S-sums and Z-sums can be 
transformed into each other, a quasi-shuffle algebra can be carried over to S'-sums. 

We remark that the quasi-shuffle algebra can be derived directly (without using the notion of 
Z-sums) as follows. Consider an alphabet A, where pairs (m, x) with m £ N and x 6 M* form 
the letters, i.e., we identify a S-sum 



{xi,x 2 ,. ■ -,x k ;n) 



with the word (mi,xi)(m2,X2) ■ ■ ■ (mfc, Xfc)- We define the degree of a letter (m,x) £ A as 
| (m, x)\ = m and we order the letters for mi, m 2 £ N and xi, X2 € M* by 



(mi,xi) 
(mx,xi) 
(mi, -xi' 



-< (m 2 ,x 2 ) 
■< (mi,x 2 ) 
~< (mi,xi) 



if 
if 
if 



mi < 777-2 
|xi| < |x 2 | 

xi > 0. 



We extend this order lexicographically to words. Using this order, it can be shown analogously 
as in [5j (compare |71j ) that the S'-sums form a quasi shuffle algebra which is the free polynomial 
algebra on the Lyndon words with alphabet A, see also |104| . 

As a consequence the number of algebraic independent sums in this algebra can be counted 
by counting the number of Lyndon words. If we consider for example an alphabet with n letters 
and we look for the number of basis sums with depth d, we can use the first Witt formula [135\ 
11361 HU6| : The number of Lyndon words of length d over an alphabet of length n is given by 



N A (d) 



q\d 



n' 



where 

'1 if n = 1 

n(n) = < if n is divided by the square of a prime 

(— l) s if n is the product of s different primes 

is the Mobius function [93, 69J. In the following we call the relations from the quasi-shuffle algebra 
also stuffle relations; cf. |28j. 

In the subsequent example we look at S-sums of depth 2 on alphabets of 4 letters. Hence 
we obtain | J2q\2 f 1 (§) 4 9 = 6 basis sums. We can use an analogous method of the method 
presented in [2H H] for harmonic sums to find the basis S-sums together with the relations for 
the dependent S-sums. 

For example, consider the letters (1, ^), (1, — |), (3, |), (3, —\)- At depth d = 2 we obtain 16 
sums with these letters. Using the relations 



Sl,3 



,1 1 

.3 V2' 2 

Sl,3 (~2J ~\ 
Sl,3 (21 ~\ 
S ~< 11 



1 

2' 



'1,3 V~2' 2 

Si,i( 

s 

Si,i 

s 3)3 (- _ 

^3,3 

S3,: 



1,1 U> 2 
1 _1 
2' 2 



1 

2' 



|;n 



2 ' 2 ' 71 



'i _i. 

v2' 2' 



7? 



-S 3 ,i 
-S 3 ,i 
-S 3 ,i 
-S 3 ,i 

|St(- 

IS1 (2! 



l ,i;n) +Si Q;n) S 3 (|;n) + S 4 Q;n) 



\ ,-|;n)+S 



2' 



77 



S 3 + S 4 



-|, |;n) +Si (|;n) S 3 (-|;n) +S4 
3,-i;n) + Si (-5577) S 3 (i;ra) + S 4 (-|;«) 

5 ;n) 2 + 2 S 2 (|;n) 



77 



+ |S 2 (i;n) 



-Si,i (-5, 3; 77) + Si (-|; n) Si (§; 77) + S 2 (-|; n) 



1. 



1 ;n) 2 + 5 s 6(i;n) 



2 

5S3 (2; ra ) 2 + sS 6 Q;n) 



-S 3 , 3 (-5, 5; n) + S 3 (-|; n) S 3 Q; n) + S 6 (-|; 77) 
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we find the 6 basis sums 

S3,i ibb n ) ' S 3,i {-hh n ) ' S 3,i ib~b n ) ' S 3,i {-b~b n ) ' 

Si,i (-2> b n ) ,s 3j 3 (- 2 , b n ) > 

in which all the other 10 sums of depth 2 can be expressed. 
9.2 Differential Relations 

In Section|8]we described the differentiation of S-sums with respect to the upper summation limit. 
Similarly to the application of differentiation to harmonic sums, cf. 27, 2(i|. and cyclotomic sums, 
cf. [5], the differentiation leads to new relations in the class of S-sums. For instance we find 

^-S 2 (2; n) = -S 3 (2; n) + H (2) S 2 (2; n) + H , .-i(l) + 2H 0)0 ,i(l) + H ,i,_i(l) . 
an 

Subsequently, we collect the derivatives with respect to n : 

S ai ,...,a k (bi, h;n)W = |^S 0l ,..., Bh (6i, . . . , b k ;n) ;N € n| . (9.2) 
Continuing the example from above we get the additional relations: 

s 3 ,i(i =^^Mb n ) - \§^ 2 ^bb n ) - l R ho{l)Mb n ) 

+ Ho(3) S 2 ,i(b \;n) - 2~ H i(i) H o,i,o(2) + ^ H o,o,i,o(i) 

s 3 ,i(-±,-§;n) = i^^ S3 (i ;n ) " \§^ 2 ^~b~b n ) + ^ H _i, (i) S2 (-5; n ) 

+ 2 H °(i) ^,i{-\,-\;n) - 2 H _i(i) H -i, (^) 

" ^ H -i A i,o(l) + ^ H OAi,o(l) " ^H_i(|) H (i) S 2 (-|;n) 

^,i{-bb n ) =\^z{-b n ) " \^ s ^(-bb n ) - \^ s 2,i{b-b n ) 

-\*lfl{l) S 2(-2^)+^H_i Q)S 2 (i;n) 

+ H (|) S 2jl (-i, i; n) + ^Ho(l) S 2) i(|, -|;n) 

- S 3 ,i(i -|;n) + ^Hi (|) Ho,-i,o(i) + (1) Ho,i,o(£) 

+ ^_i A _ 10 (i) - ^H OA -i,o(|) - ^Ho(i) Hs(J) S 2 (I;n) 

" 2 H io,-i,o(i) " 6 Ho ^ Ss (-2 ; n ) " 2 S2 ^ ! n ) S2 (^' n ) • 

Using all the relations we can reduce the number of basis sums at depth d = 2 to 3 by introducing 
differentiation. In our case the basis sums are: 

S3,l[b~b n ) > S i,i( _ I> b n ) > S 3,3(-§, b n ) ■ ( 9 - 3 ) 

Note that we introduced the letters (2, ±|) and (3, ±j), however these letters appear just in sums 
of depth 2 with weight 3 and are only used to express sums of weight 4. 
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9.3 Duplication Relations 

As for harmonic sums and cyclotomic sums, cf. |1231 l2T l l26 ll28l l5]. we have a duplication relation: 
Theorem 9.1 (Duplication Relation). For a% G N, b{ G M + and n G N we have 

^ S 0mj ... jai (±6 TO , . . . ,±&i;2 n) = Qi _ m Sa m ,...,ai (&m> • • • 

where we sum on the left hand side over the 2 m possible combinations concerning ±. 
Proof. We proceed by induction on m. For m = 1 we get 

2d 



S (-6,2n) = T { -^ = Y (-^---^—] =± Sa (b 2 ;n)-Y 7 ± 



2i - l) a 



= l Sa (6 2 ;n)-V ( 7 V + / ^-tV) 

2 a v ; ^V( 2 *) a (2i — l) a (2i)V 

= i-S a (6 2 ;n)-S a (6;2n) + ^S a (6 2 ;n) 
= ^lS a (6 2 ,n) -S a (6;2n). 
Supposing that the theorem holds for m > 1, we obtain 
^ S 0ro+ i,...,ai(±&m+l J • • • , ±6i; 2n) = 

=E(^S^Es«(±^,..,±6 i; 2i) 

±6 2i_1 \ 

+ (2i _ T)t +1 E s<w-. ai (±^ • • • , ±h; 2 i + 1) j 

= 2 E J^Ll E Sa m ,..., ai (±& m , . ... ±615 2 i) 

= 2 E (2i)«Li 2 E^iai-"» Sam -- ai ' ' • ' ') 

i=l 

^^r~7~~77Sa m+ i,a m ,...,ai(i&m+l) i&mj ■ • • > i&lj n ) 



□ 

Considering again the example from above, the duplication relations are 

S 3 , 3 (|, b 2n) =^S 3 , 3 (i |; n) - S 3 , 3 (^, ^; 2n) - S 3 , 3 (^, i; 2n) - S 3)3 (±, ^; 2n) 

Si,i(|,|;2n) =S lil (i,i;n)-Si, 1 (^,^;2n)-Si,i(^,i;2n)-Si,i(i ^;2n) 
S 3 ,i(|,i;2n) =|S 3 ,iQ,|;n) - S 3 ,i(^, ^; 2n) - S 3)1 (^, ^; 2n) - S 3 ,i(i, ^; 2n) 
Si )3 (|, |; 2n) =|Si, 3 (|, |; n) - S lj3 (^, ^;2n) - Si, 3 (^, §; 2n) - S lj3 (±, ^; 2n) . 

Note that from duplication we do not get a further reduction in this particular example. We 



would have to introduce new sums of the same depth and weight to express the basis sums (9.3) 
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Index set 


sums 


basis sums 


dependent sums 


{ai,ai},{xi,xi} 


1 





1 


{a 1 ,a 1 },{xi,x 2 } 


2 


1 


1 


{oi,o 2 },{x 1 ,a? 1 } 


2 


1 


1 


{a 1 ,a 2 },{x 1 ,x 2 } 


4 


2 


2 



Table 1: Number of basis sums for different index sets at depth 2. 



Index set 


sums 


basis sums 


dependent sums 


{a\ 


ax,ax},{xx,xx 


Xx} 


1 





1 


{a x 


ax,ax},{xx,xx 


x 2 } 


3 


1 


2 


{ax 


ax,ax},{xx,x 2 


X3} 


6 


2 


4 


{ax 


ax,a 2 },{xx,xx 


Xx} 


3 


1 


2 


{ax 


ax,a 2 },{xx,xx 


x 2 } 


9 


3 


6 


{ax 


ax,a 2 },{xx,x 2 




18 


6 


12 


{ax 


a 2 ,a 3 }, {xx,xx 


Xx} 


6 


2 


4 


{ax 


a 2 ,a 3 },{xx,xx 


x 2 } 


18 


6 


12 


{ax 


(12,(13}, {xx,x 2 


X3} 


36 


12 


24 



Table 2: Number of basis sums for different index sets at depth 3. 



9.4 Examples for Specific Index Sets 

In this Subsection we want to present the numbers of basis sums for specific index sets at special 
depths or weig hts. In the Tables [TJ [2] and |] we summarize the number of algebraic basis sums 
at the possible index sets at depths 2, 3, 4 respectively. To illustrate how these tables are to be 
understood we take a closer look at the depth d = 3 example with the index sets {ax, ax, 0,2} and 
{xx, x 2 , X3)}; these sets stand for multi-sets and each element is allowed to be taken once to build 
all the possible S'-sums. There we obtain the 18 sums: 

(x 1 ,x 2 ,x 3 ;n) , 

u ai ,ai ,a2 u ai ,ai ,0,2 

(x 2 ,x 3 ,x 1 ;n) , (x 3 ,x 1 ,x 2 ;n) , (x 3 ,x 2 ,x 1 ;n) , 

u ai .a>2 ,a\ (x 1 ,X2,x 3 ;n) ,S 

(L\ .Cl2;Q>l 

(x!,x 3 ,x 2 ;n) ,S (ll ,Cl2 ,Q>1 

(x 2 ,x 3 ,xi,n) , 

u ctl ,02, 01 

(x 3 ,xx,x 2 ;n) , (x 3 ,x 2 ,x 1 ;n) , 

Sa2,ai,ai(*^l: 3*2 , X 3 ; TlJ ; S a2i a liai \ X\ , X 3: X 2l TlJ , S 02 itll itll [x 2 , X\, X 3 , TlJ , 
Sa2,ai,ai(*^2: Xx, TlJ , S a2i a liai (x 3 , X\, X 2l TlJ , S 02 itll itll (^3 , X 2 , Xx, TlJ . 

Using all applicable algebraic relations (provided by the quasi-shufne algebra) like, e.g., 

S Q1)Q1) a 2 (iEi,,:E2, x 3l n) = S a2 (iE3; n) S ai , ai (xi, x 2l n) + S auai +a 2 (xi,x 2 x 3 ;n) 

-S ai (xr,ri) Sa. 2 , ai {x 3 ,x 2 ;n) - S a2]2ai (x 3 ,xix 2 ;n) + S a2 ^ auai (x 3 ,x 2 ,xi;n) 

we can express all these 18 sums using the following 6 basis sums (with the price of introducing 
additional S-sums of lower depth subject to the relations given by the quasi shuffle algebra): 

^(12 ,ai ,ai 

(x 1 ,x 2 ,x 3 ;n) , 

^02 ,ai ,oi 

(x 1 ,x 3 ,x 2 ;n) , (x 2 ,x 1 ,x 3 ;n) , 

Sa2,ai,ai(*^2, X 3: Xx, TlJ , S a2i a liai (x 3 , Xx, X 2l TlJ , S 02 itll itll (^3 , X 2 , Xx, TlJ . 

In Table [4] we summarize the number of algebraic basis sums at specified weights for arbitrary 
indices in the Xi , while in Table [5] we summarize the number of algebraic and differential bases 
sums for x% 6 {1, —1, 1/2, —1/2, 2, —2} and where each of the indices {1/2, —1/2, 2, —2} is allowed 
to appear just once in each sum. To illustrate how these tables are to be understood, we take 
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Index set 


sums 


basis sums 


dependent sums 


{ai,ai,ai,ai}, {xi,xi,xi,xi} 


1 





1 


{ai,ai,ai,ai}, {xi,xi,xi,x 2 } 


4 


1 


3 


{ai,ai,ai,ai}, {xi,xi,X2, x 2 } 


6 


1 


5 


{ai,ai,ai,ai}, {xi,xi,x 2 ,x 3 } 


12 


3 


9 


{ai,ai,ai,ai}, {xi,X2,X3,X4} 


24 


6 


18 


{ai, ai, ai, 02}, {xi, xi, xi, xi} 


4 


1 


3 


{ai,ai,ai,a 2 }, {xi,xi,xi,x 2 } 


16 


4 


12 


{ai, ai, ai, 02}, {xi,xi,X2, X2} 


24 


6 


18 


{ai,ai,ai,a 2 }, {xi,xi,x 2 ,x 3 } 


48 


12 


36 


{ai, ai, ai, a 2 }, {xi,x 2 ,x 3 , x 4 } 


96 


24 


72 


{ai,ai,a 2 ,a 2 },{xi,xi,xi,xi} 


6 


1 


5 


{ai,ai,a 2 ,a 2 }, {xi,xi,xi,x 2 } 


24 


6 


18 


{ai, ai, a 2 , a 2 }, {xi,xi,x 2 , x 2 } 


36 


8 


28 


{ai,ai,a 2 ,a 2 }, {xi,xi,x 2 ,x 3 } 


72 


18 


54 


{ai, ai, a 2 , a 2 }, {xi,x 2 ,x 3 , x 4 } 


144 


36 


108 


{ai,ai,a 2 ,a 3 }, {xi,xi,xi,xi} 


12 


3 


9 


{ai, ai, a 2 , a 3 }, {xi,xi,xi,x 2 } 


48 


12 


36 


{ai, ai, a 2 , a 3 }, {xi,xi,x 2 , x 2 } 


72 


18 


54 


{ai,ai,a 2 ,a 3 }, {xi,xi,x 2 ,x 3 } 


144 


36 


108 


{ai, ai, a 2 , a 3 }, {xi,x 2 ,x 3 , x 4 } 


288 


72 


216 


{ai,a 2 ,a 3 ,a 4 }, {xi,xi,xi,xi} 


24 


6 


18 


{ai,a 2 ,a 3 ,a 4 }, {xi,xi,xi,x 2 } 


96 


24 


72 


{ai, a 2 , a 3 , a 4 }, {xi,xi,x 2 , x 2 } 


144 


36 


108 


{ai,a 2 ,a 3 ,a 4 }, {xi,xi,x 2 ,x 3 } 


288 


72 


216 


{ai, a 2 , a 3 , a 4 }, {xi, x 2 , x 3 , x 4 } 


576 


144 


432 



Table 3: Number of basis sums for different index sets at depth 4. 



a closer look at two examples. At weight w = 3 with index set {xi,X2,x 3 }, we consider the 19 
sums: 

S 3 (xix 2 x 3 ; n) , Si )2 (a;i, x 2 x 3 ; n) , Si i2 (x 2 , XiX 3 ; n) , Si ;2 (xiX 2 , x 3 ; n) , Si j2 (x 3 , XiX 2 ; n) , 
Si, 2 (xix 3 , x 2 ; n) , Si, 2 (x 2 x 3 , Xf,n), S 2 ,i(xi, x 2 x 3 ; n) , S 2; i(x 2 , xix 3 ; n) , S 2; i(xix 2 , x 3 ; n) , 
S2,i(x 3 , xix 2 ; n) , S 2 ,i(xia;3, x 2 ; n) , S 2j i(x 2 x 3 ,xi;n) , Si,i,i(xi, x 2 , x 3 ; n) , Si,i,i(xi, x 3 , x 2 ; n) , 
Si,i,i(x 2 ,xi,x 3 ;n) , Si,i,i(x 2 ,x 3 ,xi; n) , Si,i,i(x 3 , xi, x 2 ; n) ,Si t \,\{x 3 ,x 2 ,x\\n) . 

Using all applicable algebraic relations (provided by the quasi-shume algebra) like, e.g., 

Si, 2 (x 2 x 3 , xi;n) = Si(x 2 x 3 ; n) S 2 (xi; n) + S 3 (xix 2 x 3 ; n) - S 2 ,i(xi, x 2 x 3 ; n) , 
Si,i,i(xi,x 2 ,x 3 ;n) = -Si(x 3 ;n) Si,i(x 2 ,xi;n) + Si(xi;n) Si i i(x 2 ,x 3 ;n) 

+S 2 ,i(xix 2 , x 3 ; n) - S 2 a(x 2 x 3 , xi;ri) + Si ; i i i(x 3 , x 2 , Xi; n) 

we can express all these 19 sums using the following 9 basis sums: 

S 3 (xix 2 x 3 ; n) , S 2 ,i(xi, x 2 x 3 ; n) , S 2 ,i(x 2 , xix 3 ; n) , S 2 ,i(xix 2 , x 3 ; n) , S 2 ,i(x 3 , xix 2 ; n) , 
S 2 ,i(xix 3 ,x 2 ;n),S 2i i(x 2 x 3 ,xi;n),Si ; i i i(x 3 ,xi,x 2 ;n),Si ; i i i(x 3 ,x 2 ,xi;n) . 

At weight w = 2 where each of the indices from the index stet {1/2, —1/2, 2, —2} is allowed to 
appear just once, we consider the 38 sums: 

S_ 2 (n) , S 2 (n) , S_i,_i(n) , S_i,i(n) , Si,_i(n) , Si,i(n) , S 2 (-2; n) , S 2 (-|; n) , 
S 2 (i; n) , S 2 (2; n) , S M (-2, -1; n) , S M (-2, -\- n) , S M (-2, \; n) , Si,i(-2, 1; n) , 
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weight 


index set 


sums 


basis sums 


dependent sums 


2 


{xi,xi} 


2 


1 


1 




{xi,x 2 } 


3 


3 


1 


3 


{x!,Xl,Xl} 


6 


3 


3 




{x 1 ,x 1 ,x 2 } 


12 


6 


6 




{x 1 ,x 2 ,x 3 } 


19 


9 


10 


4 


{x 1 ,x 1 ,x 1 ,x 1 } 


20 


8 


12 




{x 1 ,x 1 ,x 1 ,x 2 } 


50 


20 


30 




{x 1 ,x 1 ,x 2 ,x 2 } 


64 


24 


40 




{xi,xi,x 2 ,x 3 } 


106 


40 


66 




{x 1 ,X 2 ,X 3 ,X4} 


175 


64 


111 


5 




70 


25 


45 




{xi,Xl,Xl,Xl,X 2 } 


210 


70 


140 




{x 1 ,x 1 ,x 1 ,x 2 ,x 2 } 


325 


105 


220 




{xi,xi,xi,x 2 ,x 3 } 


555 


175 


380 




{xi,xi,x 2 ,x 2 ,x 3 } 


725 


225 


500 




{xi,xi,x 2 ,x 3 ,x 4 } 


1235 


375 


860 




{x 1 ,X 2 ,X 3 ,X4,X 5 } 


2101 


625 


1476 


6 


{x 1 ,x 1 ,x 1 ,x 1 ,x 1 ,x 1 } 


252 


75 


177 




{xi,Xl,Xl,Xl,Xl,X 2 } 


882 


252 


630 




{x 1 ,x 1 ,x 1 ,x 1 ,x 2 ,x 2 } 


1596 


438 


1158 




{x 1 ,x 1 ,x 1 ,x 2 ,x 2 ,x 2 } 


1911 


522 


1389 




{x 1 ,x 1 ,x 1 ,x 1 ,x 2 ,x 3 } 


2786 


756 


2030 




{x 1 ,x 1 ,x 1 ,x 2 ,x 2 ,x 3 } 


4431 


1176 


3255 




{x 1 ,x 1 ,x 2 ,x 2 ,x 3 ,x 3 } 


5886 


1539 


4347 




{X 1 ,X 1 ,X 1 ,X 2 ,X 3 ,X4} 


7721 


2016 


5705 




{X 1 ,X 1 ,X 2 ,X 2 ,X 3 ,X4} 


10251 


2646 


7605 




{X 1 ,X 1 ,X 2 ,X 3 ,X4,X 5 } 


17841 


4536 


13305 




{xi,X 2 ,X 3 ,X4,X 5 ,X 6 } 


31031 


7776 


23255 



Table 4: Number of basis sums for different index sets up to weight 6. 



Si,i(-2,2;n),Si,i(-l,-2;n),Si,i(-l,-|;n),Si,i(-l,|;n),Si,i(-l,2;n) ) 
Si,i(-5,-2;n) ,Si,i(-|,-l;n) , Si,i(— |, \; n) , S M (-i, 1; n) , Si,i(-|, 2; n) , 
S M (|,-2;n) ,Si,i(|,-l;n) ,S M (|,-|;n) ,S M (|,l;n) ,Si,i(§,2;n) , 
Si ; i(l,-2;n),S 1;1 (l,-|;n),S 1 , 1 (l,i;n),S 1 ,i(l,2;n),Si, 1 (2,-2;n),Si, 1 (2,-l;n), 
Si,i(2,-i;n),Si il (2,i;n),S 1 ,i(2,l;n). 

Using all algebraic and differential relations like, e.g., 

Si,i(-|,2;n) =- ^S_i(n)+H_ li0 (l)+Si(-i;n) Si(2; n) - S M (2, -\- n) 
Si,i(-|,-2;n) =-^Si(n)-H li0 (l) + Si(-2;n) Si(-§; n) - S M (-2, -\; n) 
we can express all these 38 sums using the following 17 basis sums: 

Si,-i(n) , Si,i(-2, -I; n) , S M (-2, \;n) , Si,i(-2, 2; n) , S M (-1, -2; n) , 
Si,i(-l,-|;n) ,Si ; i(-l, \;n) , Si i i(-1, 2; n) , Si i i(-|, |; n) ,S\,i(\,-\;n) , 
S 1)1 (l,-2;n),S 1)1 (l,-i;n),S 1 ,i(l,i;n),S M (l,2;n),Si,i(2,-2;n), 
Si,i(2,-|;n),Si,i(2,|;n). 
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weight 


sums 


a-basis sums 


d-basis sums 


ad-basis sums 


1 


6 


6 


6 


6 


2 


38 


23 


32 


17 


3 


222 


120 


184 


97 


4 


1206 


654 


984 


543 


5 


6150 


3536 


4944 


2882 


6 


29718 


18280 


23568 


14744 



Table 5: Number of basis sums up to weight 6 with index set {1, —1, 1/2, —1/2, 2, —2}. Each of 
the indices {1/2, —1/2, 2, —2} is allowed to appear just once in each sum. 



We remark that the found relations are induced by the quasi-shuffle algebra. Using the 
summation package Sigma |113j and the underlying difference field theory [TTJ 1118} 1116} 1117] we 
could verify that the found basis sums are also algebraic independent in the ring as sequences; 
for related work see also [68]. To be more precise, the given 5-sums evaluated as sequences form 
a polynomial ring which is a subring of the ring of sequences. 

In general, we obtain the following relations 

N D (w) = N s (w) - N s (w - 1) 
N AD {w) = N A (w) - N A {w - 1) 

where Ns(w), Nd(w) and N A d{w) are the number of sums, the number algebraic basis sums and 
the number of basis sums using algebraic and differential relations at weight w respectively. 



10 Relations between 5-Sums at Infinity 

In this Section we will state several types of relations between the values of iS-sums at infinity 

Smi,...,m p (%l, ■ ■ ■ ,x p ;oo) (10.1) 

(or equivalently between the values of generalized polylogarithms) that extend various relations 
for infinite harmonic sums given in [28]. This enables one to reduce the expressions that arise 
during the calculations of, e.g., the (inverse) Mellin transform, of the differentiation of S-sums, 
or of asymptotic expansions given below. We remark that the relations of S*-sums with X{ being 
roots of unity have been considered in [5] already. 

Stuffle relations: The first type of relations originates from the algebraic relations of 5-sums, 
see Section 9.1 Those relations which contain only convergent 5-sums remain valid when 
we consider them at infinitjj^j We will refer to these relations as the quasi shuffle or stuffle 
relations. 



Duplication relations: The duplication relations from Section 9.3 remain valid if we consider 
sums which are finite at infinity, since the limits n resp. 2n — > oo agree. 

Shuffle relations: We can generalize the shuffle relations form [123J for harmonic sums to S- 
sums. For convergent sums we have: 

Smi,...,m p (xi, ■ ■ ■ , x p ] oo) S klt ...,k q {yi, ...,y p ;oo) = 

lim y\'^'m 1 ,...,m p {xi, • • • , Xp] n — i) Sfc 2 . . ^(^ • • • > Vp', i) 



n— >oo 

i=l 



6 In general, relations containing sums of logarithmic growth such as Sx(n) can be also formally utilized in 
asymptotic expansions. However, in the following only convergent sums are considered. 
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Using 

lim \^ y\^rni,...,m p ( x li ■ ■ ■ > x p'i n — i) Sk 2! .., 1 k q (y2i • • • ; Up] i) 

n— s>oo ' i^l 
i=l 

E l fk± + mi — 1 — a\ x| 
1 mi - 1 / ^ jroi+fci-a 

a=l v 1 7 i=l 



A S m2i ..., mp (x 2 , ...,x p ;i-j) S k2t ... )kq (y 2 , y P ;j) 



i=i J 

m i /j , i \ n i 

l mi - 1 / ^ ^mi+fci-o 

a=l x 7 i=l 



* (fr) ^ $k 2 ,...,k q {y2, ...,y p ;i-j) S m2i ..., mp (x 2 , . . . ,x p ; j) 



3=1 J 

we can rewrite the right hand side in terms of S'-sums. We will refer to these relations as 
the shuffle relations since one could also obtain them from the shuffle algebra of generalized 
poly logarithms . 

Duality relations: We can use duality relations of generalized polylogarithms; compare [28] 
where these relations are considered for harmonic sums. 

1 — x — > x: For a generalized polylogarithm H 01(02( ... )0fc _ 1)afc (l) with a& ^ this leads to the 
relation: 

H 0l ,oa,...,ofe_i,ofc(l) = Hi_ afc) i_ afe _ lr .. ! i_ a2! i_ ai (l) . 
E.g., we get H2(l) = H_i(l) which gives in S-sum notation the new relation 

Si(|;oo) = -S_i(oo) . 
— > x: This transform leads, e.g., to the relation rWl) = — H_i(l) + H_i (f) which gives 

L+X 3 

in S-sum notation again the relation 

Si(i;oo) = -S_i(oo) . 
— > x: For c, d G K, d / —1 this transform can be calculated similarly as described in 



3.3 



The transform — > x produces, e.g., the relation H_2(l) = |(H_s(l) + 



Section 

H i (1)) which gives in S-sum notation the new relation 

Si(-|;oo) = -S_i(oo) + Si(-^;oo) . 



If one considers a certain class of infinite S'-sums (10. 1) (with finite values) with given weight 
and a finite index set from which the x% can be chosen, one can try to determine all relations 
from above that stay within the specified class. In general, we first have to extend the alphabet 
generated by the relations used. Then one exploits all relevant relations that lead back (after 
several applications) to the original alphabet. 

In the following Subsections we apply this tactic to a series of alphabets. We summarize in 
Table [3] how all the so far calculated relations for the particular alphabets contribute in reducing 
the number of infinite S-sums to a smaller set of sums. For the found constants we searched 
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alphabet 


weight 


finite sums 


dependent sums 


basis sums 




1 


1 





1 


{1,-1} 


2 


4 


3 


1 




3 


12 


11 


1 




4 


36 


35 


1 




1 


1 





1 


{!>*} 


2 


4 


2 


2 




3 


12 


8 


4 




4 


36 


24 


12 




1 


1 





1 


{1,1,2} 


2 


5 


3 


2 




3 


18 


14 


4 




4 


66 


53 


13 




1 


2 


1 


1 


{1,-1, 1} 


2 


9 


7 


2 




3 


36 


29 


7 




4 


144 


113 


31 




1 


3 


1 


2 


S\ _i I _I\ 

X i ■*•> 2' 2-1 


2 


16 


12 


4 




3 


80 


63 


17 




4 


400 


301 


99 




1 


3 


1 


2 


{1, —1, 2 ' ~~ 2 ' 2, —2} 


2 


20 


16 


4 




3 


124 


107 


17 




4 


788 


678 


110 



Table 6: Number of basis constants (convergent infinite 5-sums) for different alphabets up to 
weight 4. 

for new relations using PSLQ [59J using representations of up to 2000 digits. No further relations 
were found. These constants, which we also call basis constants, enlarge the number of constants 
generated by the multiple zeta values [28 . They are further extended by other special numbers 
related to cyclotomic harmonic sums and their generalization at infinity, cf. e.g. |37[ \5\ . 

Subsequently, the found basis constants for the concrete alphabets are given. In the simpler 
cases also the relations between the dependent sums are printed explicitly Q Given an expres- 
sion with the corresponding weight and alphabet and given these relations, one can bring the 
expression to its reduced form, i.e., in terms of the basis constants. Here one needs to replace 
the S-sum on the left hand side of a given relation by its right hand side. 

10.1 Alphabet {1, |} : 

w=l: 

Basis constant: S\ (|;oo) = log(2) 
Relation: — 

w=2: 

Basis constants: 52 (oo) = £2, Si t i (5, 5, °°) 

7 In http://www.risc.jku.at/research/combinat/software/HarmonicSums/ all relations are available up to 
weight four. The most involved calculation is the alphabet Xi £ {1, — 1, 1/2, —1/2, 2, —2}; the relations are obtained 
executing ComputeSSumlnf Basis [4, {1, — 1, 1/2, —1/2, 2, — 2}, ExtendAlphabet — > True] in about 36 hours. 
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Relations: 



<Si,i ( 2' 1;o ° 



s 2 I ^ 



2 



Basis constants: S 3 (oo) = C3, £1,2 (5, 5; 00) , 51,1,1 (5, §, §; 00) , Si, 1,1 (|, 
Relations: 

log 3 (2) , 1, ,_ v 7 Cs^o Q ^ 
= — 7. 1- 7, !og(2)C2 - ~pr + 351,2 



^3 I 2^ 



£2,1 ( 1, 2-; 00 

S 2 ,i(oo 



5*1,2 (^2' 1;o ° 
Si,i,iQ,l^;oc 



51,1,1 ( 2' x ; 00 



2 Alphabet {1, |, 2} : 



6 



--(log(2)C 2 ) + C3 
^ + ^o g (2)C2 + | 



2C3 

5C3 



3C3 

4 



Basis constant: S± (^; 00) = log(2) 
Relation: — 



2: 

Basis constants: ^2(00) = C2, Si,i (|, |; 00) 
Relations: 



<H - -2 iog2 ^ + l 

5i,iQ,i;oo) = I 

5i,iQ,2;oo) = ^ 

3: 

Basis constants: S 3 (oo) = C3, <Si,2 (5, \; 00) , Si, 1,1 (5, 5, |; 00) , Si, 1,1 (\, 
Relations: 

2 ;o ° = -2 C2log(2) + -e— + T" 
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5 2 ,i ( 2' 1;o ° 



s 2 ,i ( 2' 2 ; 00 



5 2 ,i ( 1, ^; 00 
52,1(00 

5-i 2 ( 2" hoc 



5i, 2 l^ 2 ! 00 
5l,l,l ( -, 1, 1; 00 



5i,i,i y-, 1,2; 00 

/ll 

Sx,!,! ^2> 2 '2 ;0 ° 

51,1,1 [2' 2 ' 1; 00 



1 1 

2' 2 



ic 2 log(2) + ^^ + 35 li2 
-^C2log(2) + C 3 

i ( _ 3)&log(2) + ?lM 
i (2log(2) + !^ + | 

2C3 

5C(3) 
8 

^C 2 log(2) 

-ic 2 log(2) + ^^ + log(2)5 1 , 1 



; 00 



7Cs 



1 1 

2'2 ; 



2S U 



1 



= -C 2 log(2) + 5 li2 

4 

7C(3) 



|c 2 Iog(2)-| 
fc 2 log(2) + ^ 



1 1 

2'2 5 



00 — 



5Cs 
12 



11 V 7( 3 

-, -;oo H 



2' 2' 



24 



10.3 Alphabet {1, —1, |} : 



w=l: 



Basis constant: 5_i(oo) = — log(2) 
Relation: 



-5_i (00) 



w=2: 



Basis constants: 5 2 (oo) = ( 2 , 5i ; i (5,5; 00) 
Relations: 



5- 2 (00) 
5 2 Q;oo 
5-i -1(00) 



5j.. ( -i^; 00 



5-i, 1(00) 



2 

C 2 log 2 (2) 
2 2 

C2 log 2 (2) 
2 2 



C 2 log 2 (2) 
2 



2 +2 M^ ;o ° 



C2 log 2 (2) 
2 2 
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/ 


1 






4 


2 




oo 






l 






( 




oo 






1 




( 


2' 2; 


oo 



2 

3^ 
2 



w=3: 

Basis constants: 5 3 (oo) = (3, 5i,2 (-1, \\ 00) , Si,i,i (-1, |, \\ 00) , Si,i,i (-1, \, 1; 00) 

5l,l,l (2, -1, \\ 00) , S*!,!,! (5, -1; 00) , S*!,!,! ( 2 , 2l \\ °°) 

Relations: 



5-3(oo) = -1(3C 3 ) 



(I \ log(2) 3 log(2)C 2 7C 3 



5'-2-i(oo) = 



31og(2)C 2 5C 3 
"28 

o ( 1 1 "\ log(2) 3 log(2)C2 7C 3 1 / 1 \ 

52 H" '2'°°J = ^8 3~ + 72 + 3*M~ '2' °° J 



5- 2 ,i (00) = -1(5C 3 ) 



, 8 

/l \ 41og(2) 3 log(2)C 2 4C 3 2 /1 \ 

5 2 ,i^,-l;ooj = — — -_-- 5ll2 ^ -l,-;ooj 

5 2 ,i Q,*; 00) = -I(log(2)C 2 ) + C3 + 5 li2 (-l,i;oo) 

5 2 ,iQ,l;oo) = -l(log(2)C 2 ) + C3 

$2,-1(00) = -I(31og(2)C 2 ) + | 

/ 1 \ log(2) 3 log(2)C 2 C3 

S 2 ,i(oo) = 2C 3 

5-1,-2(00) = -log(2)C 2 + 1:K ' ! 



8 

5-1,2(00) = Cs 



log(2)C 2 



5i, 2 Q,-l;oo) = -^(13C 3 ) 

log(2) 3 log(2)C 2 . 13(3 , l q ( 1 1 \ 



1 1 

5i, 2 ( 2' 1;0 ° 



log(2) 3 log(2)C2 Cs 
>-i,-i,-i(oo) = 



o f 1 1 1 \ 171og(2) 3 131og(2)C 2 6IC3 , , Q fl 1 

5i,i,i^-l,-l,-; 00 j = ~^^ + g ^- 41 ° g ^U'2 ;0 ° 

-^5i, 2 
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log(2) 3 log(2)C 2 7C 3 
5_i_i,i(oo) = + — 



6 2 4 

1 1 



-l;oo) = -^(21og(2) 3 )-31og(2)C 2 + ^+41og(2)5 1 , 1 ( 



i,i ( 2'2 ; °° 



+45io ( -l,x;oo 



2' 

c , x log(2) 3 , log(2)C2 , Cs 
_i(oo) = g — + g + ~~ 



D Z O 

c f , , !. J 51og(2) 3 51og(2)C 2 6IC3 4 / 1 \ 
5i,i,i^-l,l, 2 ,ooj - — - + _ + -5 1>2 ^-l,-,ooJ 

log(2) 3 log(2)C 2 7C 3 
S-i,i,i(oo) = ^ + ^7~ "Y" 

c f 1 i -1 . "\ 71og(2) 3 41og(2)C 2 35C 3 . , 9W (I 1 

-^i, 2 (-i,^;oo) 

c f 1 11.^ 21og(2) 3 log(2)C 2 155(3 1 ( 1 \ 
Si,i,i^ 2 ,-l,l,ooj - ^^^ + — ____5 li2 ^_ l, 2 ,ooj 



-Jsin ( -l,-,l;oo 
2 ' ' I 2 



(I 1 \ log(2) 3 71og(2)C 2 217C 3 . WoW /l 1 
Si,i,i^ 2 , 2 ,l,ooj = — ^ + l4T +1 ° g(2)5l ' 1 U'2' ( 

+ ^1,2 ^-1, ^ °o^ + isi >M ^-1, ^, 1; oo^ 



5 M ,i Q,i,^;oo 

Q f 1 1 1 ^ 3 ^ 3 

5i,i,i I -,1,1; oo J = — 



10.4 Alphabet {1, -1, \, -§} : 



w=l: 



Basis constant: 5_i(oo) = — log(2),Si (— |;oo) 



w=2: 

Basis constants: S^oo) = ( 2 , S^i (— \, \\ oo) , S^i (|, — ^; oo) , Si ; i ^; oo) 

w=3: 

Basis constants: S 3 (oo) = C3, S 3 (-^; oo^j , Si )2 f-^, oo^ , 5i )2 f-^, ^; oo , 
1 1 \ „/ ll \ „ / 1 1 \ „ / 11 



Si,i,i (-i4'-^ ;o °) H' -1 '"^ 00 ) (-^'-^ ;o °) ' 
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S.u,. ( 2'-2'-2 ; °°) ' 51,1,1 2 ; °°) ' 5w U'4' 1; °° 



10.5 Alphabet {1, -1, §, 2, -2} : 



w=l: 

Basis constant: S_i(oo) = — log(2),Si (— ^;oo) 



w=2: 

Basis constants: S 2 (oo) = C2 J £l,l Q, t,;c»^ ,<Sl,l 7,;°°^ j^l.l 

w=3: 

Basis constants: 53(00) = (3,83 ^— -;oo^ , Si,2 ^— -,— 2; 00^ , Si t 2 
11\„/11\„ / 11 







(4 










1 1 


"2'2 ; 


00 ) , 





Si,2 ( v 2'~2 ;0 ° J ' 5,2,1 V 2' 2' °° J ,Sl ' 1 ' 1 V 1 '~2'~2' 00 J ' 5,1,1,1 ' _1 '~o' o5°° 

51,1,1 ' 51,1,1 {~i~ h ~l ;o °) ' 5w H'~^ ; 

51,1,1 ("2' 2'-2 ; °°) ' 51,1,1 ("2' 2' 2 ; °°) ' 51,1,1 ("2' 2' 1; °°) ' 

51,1,1 U'-^'-^ 00 ) ' 51,1,1 U*-^ ;o °) ' 51,1,1 U'"5' 1;o °)- 



11 Asymptotic Expansion of S'-Sums 



In this Section we seek expansions^] of S-sums generalizing ideas for harmonic sums given in 
We say that a function / : R — )■ R is expanded in an asymptotic series |132j 

00 

n=l 

where a n are constants from R, if for all iV > 

iV 



n=0 V y 



x — > 00; 



G(z) 



0. When 



note that for a function g(x) we use the notion g(x) G o(G(x)) 4=> lim 

only the first coefficients of such an expansion are given (calculated), we often write f(x) ~ 
ao + a\x + a2X 2 + • ■ • + a n x n and cut off the remainder term. 

We remark that the subsequent Sections provide algorithms to expand S'-Sums which come 



close to the class of S'-Sums (for a definition see (7.6)). For the general class of S-sums the 



expansion problem is not handled completely. However, for current calculations in QCD (see, 



e.g., Section 12) the presented methods are completely sufficient. 



For further computer algebra aspects concerning the calculation of expansions we refer, e.g., to [134, 108, 701161] 
and references therein. 
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11.1 Asymptotic Expansions of Harmonic Sums Si(c; n) with c > 1 
Version 1: The asymptotic representation of the harmonic numbers Si(n) is given by [9 

S,(n) ~ 7 + login) + ^ " E ■ ("-I) 

fe=l 

where 7 is the Euler-Mascheroni constant and B n are the Bernoulli numbers |17l 11071 197] , 

fc=o ^ ' 

To compute the asymptotic expansion of Si(c;n) with c > 1 one may represent it using the 
Euler-Maclaurin formula 



Si(c;n) = 5i(c; n) + Si(ra) 



where 



Si(c;n) := ^ 



c 4 - 1 



log(c) + 



1 



n + 1 



+Ei(log(c)(n + l))+^ 



7 - log(log(c)) - log(n + 1) 

n+l 



B 



2.; 



Here Ei(z) is the exponential integral [55| [95] . 

rx c t _ 1 



+ R k (c,n) . 







dt = Ei(xlog(c)) — 7 — log(x) — log(log(c)), 



and Rk(c,n) stand for the rest term. The derivatives are given in closed form by 



d? / (*-! 
dP \ i 



d? (c l -\ 
lim — — 

i-+0 dP 



c^ iog J ' fc (c) ( 1)fc r(j + i) 



1 « — * % 

k=0 



T(j + 1 - k) 



r(j + 1) 



p 



i+l 



logj +1 (c 

i + i 



The asymptotic representation for Si (c;n) with c > 1 can now be derived using 



Ei(i) ~ exp(t) J2 



k=l 



T(k) 
t k ' 



Version 2: A more suitable approach to compute the expansion of Si(c;n) with c > 1, and 
which is implemented in the HarmonicSums package, is as follows. Differentiating the integral 
representation 

Si (c;n) 



n _J_ 

E + 

i=l 



1 _ JL 



x 



with respect to n we get 

d Si(c;n) 
dn c n 



1 x-Ho(x) + ^Hq(c) 



dx 
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Hence 



/■'x"(Ho W + Ho(c) )<fa _ Ho(e) f^ dx 
Jo x - ~ Jo x - ± 



l x n (R (x) + R (c)) d Sx(c;n) n S t (c;n) 

= oIT ^ +H (c) - . (11.2) 



/o x — - dn c n c 

We can now compute the asymptotic expansion of the integral using integration by parts. Let 



r 1 x n 

Jo 



x _i 2^ n i 



with aj 6 1 be this expansion and let Ylt=o ^ with h G M be the asymptotic expansion of Sl ^f^ . 
We note that the differentiation of this expansion yields Ylt^=2 i- e - 5 

d Si(c;n) ^ ^ -(j - 

i=2 



Plugging these expansions into equation (11.2) and comparing coefficients yields the recurrence 
relation 

b n = T^T + l^T&n-i forn > 2 
H (c) H (c) 

with the initial values bo = and b\ = j^jy- Hence we get for example 

„ i * „_ / 14174522 1091670 94586 9366 1082 150 26 6 2 2 
Si2;n ~ 2™ ^— + + _^ + + ++ + + + _ 



n iU rr n" n' n" nr nr n u n* n 



( 566733 17295 38001 1491 273 15 33 3 3 3 
l( ' n) ~ \~Jn^ + + ~16^ + 4^ + 4^ + ^ + 8^ + 2rJ + 4^2 + fri 

11.2 Computation of Asymptotic Expansions by repeated partial Integration 

Next, we deal with generalized polylogarithms of the form 

M / Hmi, W3 ,.-,m > (g) \ andM f Hbl ' b2 -- fefc(1 ~ X) ,n^ (11.3) 
\ c±x ) \ c±x ) 

with | c| > 1 where rrii G R \ (0, 1] or b { G E \ (-1,0) U (0, 1)) with / using repeated 
integration by parts. The first two lemmas show that the arguments of the Mellin transform 
under consideration are analytic. Thus iterated integration by parts is possible for c > 1 using 
the already known integral representation. The third lemma gives an integral representation for 
the special case c = 1 such that the method of repeated partial integration is also possible. The 
proofs are omitted here. 

Lemma 11.1. Let a mit m 2 ,...,m k (x) be a generalized poly logarithm with rrii G M \ (0,1] for 
1 < i < k. Then 

tt / \ H-rn 1 ,m 2 ,...,in k {x) , Hm 1 ,m 2 ,...,m k {x) — H mi;?7l2v , i|Tnfc (l) 

- n mi,m2,...,rrti- \X J , ana 

c + x c — X 

are analytic for |c| > 1 and x G (0, 1]. 
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Lemma 11.2. LetJl 

for 1 < i < k, and m& ^ 0. Then 



mi,m2,...,m(. 



(x) be a generalized poly logarithm with m,i G K \ ((—1,0) U (0,1)) 



H 



mi ,ni2 



,...,m fc (l X) 



and 



H 



mi,m2,...,m/ c 



(1-x) 



c — X 



C + X 

are analytic for |c| > 1 and x G (0, 1]. 

For c > 1, we obtain a suitable integral representation for partial integration. E.g., we get 



x + 2 ' 



x + 2 



x n x ■ H_ 2 (x) 



n x + 2 



1 l g -H_ 2 (x) 
x + 2 

1 x" 2H_ 2 (x) +x 



o n (x + 2) 2 



dx 



=H_ 3 (l)3-i/ 
n n Jo 



3 1 Z" 1 „_ 1 x(2H_ 2 (x) +x) 



(x + 2) 2 



dx 



=H- 2 (1) 
=H- 2 (1) 



3 


I 


' x n x{ 


n 




i n 






2 \ 


(" 

V n 


9n 2 J 



1 x"6x-2(x-2)H_ 2 (x) 



H- a (l) 



n (x + 2) 3 

I 1 / ,1 ^„_ 1 x( 6x-2(x-2)H_ 2 (x)) 

(x + 2)3 
1 2 



rfx 



rfx 



20 / J_ 

n 9n 2 ' 27n 3 ' 27n 4 / 9n 2 ' 9n 3 81n 4 V n 5 



For the special case c = 1 one can use the integral representation of the following lemma to 
calculate the expansion with partial integration as demonstrated above. 

Lemma 11.3. Let H m (x) = H mij?n2) ... jmfe (x) and H&(x) = H^^...^^) 6e a generalized polylog- 
arithm with rrn G R \ (0, 1] /or 1 < i < fe and h € M \ ((-1,0) U (0, 1)) for 1 < i < I 
where hi ^ 0. TTien we Ziowe 



1 — x 



1 „.n 



x n (H m (g) - H m (l)) 
1 — x 



dx 



H m (x) -H m (r 



1 — x 



dx - Si(ra)H m (l) , 



and 



where 



are finite constants. 



1 - x 

1 H m (x) - H m (l 
1-x 



1 x n H 6 (l - x) 



1 — x 



dx - H ,b(l) 



■dx, H m (l) and H ,&(1) 



Summarizing, we are able to calculate the asymptotic expansion of (11.3) using the method of 
repeated integration by parts. 



11.3 Computation of Asymptotic Expansions of S-Sums 

Subsequently, we derive an algorithm that calculates the asymptotic expansion of S'-sums 
Sai,...,a fc (&i, ■ • • , bk', n) with hi G [—1,1] and hi ^ using the technologies of the previous Sub- 
sections. 



Since these sums fall into the class of S'-sums, it follows by Section 7.2 that they can be 
represented using Mellin transforms of the form M f H "i."^W 3! ) ; n \ with m . £ R\((_i, 0)U(0, 1)) 
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and |c| > 1. Conversely, given such a Mellin transform, it can be expressed in terms of S-sums 
of the above type with bi G [—1, 1] and bi ^ 0. 

In addition, it is straightforward to see that this subclass is closed under the quasi shuffle 
product (2.2), i.e., the product can be expressed again as a linear combination of S-sums of this 
subclass over the rational numbers. 

We are now ready to present an algorithm which computes asymptotic expansions of S-sums 
S ai ,a 2 ,...,a k (bi, b2, ■ ■ ■ , b k ; n) with bi G [—1, 1] and bi / 0. The algorithm tries to attack the sums 
directly using the previous Sections. If it fails one uses the transformation x — > 1 — x from 
Section 3.2 twice and reduces the problems to S-sums that are less complicated than the input 
sum. Recursive application will finally produce the expansion of the original S-sum. To be more 
precise, the algorithm can be described as follows. 

• If S ai|( j 2) ... !afc (6i, 62, • • • , b k ; n) has trailing ones, i.e., a k = bk = 1, we first extract them such 
that we end up in a univariate polynomial in Si(n) with coefficients in the S'-sums without 
trailing ones. Expand the powers of Si(n) using ( |11.1[ ). Now apply the following items to 
each of the S'-sums without trailing ones: 



Let S 01j02j ... j0fc (6i,6 2 ,...,fcfc;n) with a k ^ 1 ^ b k , and compute mi '"^' 
the inverse Mellin transform the most complicated S-sum is S 0lja2) ... iafc (&i, 62 
express S ai ,a 2 ,...,a fe (&i, h, ■■■,b k ;n) as 



(x) 



$ai,a 2 ,...,a k (b 1 ,b2, ...,b k ;n) = M 



H 



mi,m2,...,mi 



(X) 



c + sx 



, n 



+ T 



such that in 
.,b k ;n), i.e., 



(11.4) 



where T is an expression in S-sums S^ ^ y ; (b[, 6' 2 , . . . , b' k ,; nj with b\ G [—1, 1] and b[ ^ 
(which are less complicated than S aiia2) ... )afe (&i, &2j ■ ■ • > b k ; n)) and constants. Note that 
s = ±1, c > 1 and m,i G K \ ((—1, 0) U (0, 1)) and that mi ^ 1 (see the beginning of this 
Subsection and Lemma 7.5). 



We proceed by expanding M 



-"-m^ ,m2 , • ■ ■ 

c-\-sx 



lx) 



n) in (11.4): 



— 

not all rrii ^ 1: 



m l ' m 2 1 ■ ■ ■■< rrL t 



(x) 



c+sx 



, n^j directly; see Section 



11.2 



— Transform x — > 1 — x in H m (x) as described in Section 3.2 and expand all products, 
i.e., we obtain 



\c + sx 



p 

E 

i=l 



dM 



Hfc(l 



x 



c + sx 



n + d with d, d% G 



(11.5) 



Note that due to Remark 3.5 the components of the bi are in M. \ (0, 1). 



— For each Mellin transform M 

• see Section 



-fc — - , n ) do 



!>., " : Expand f Q r+s r 
b-j = : As described in Section 



3.2 



11.2 



transform back 1 — x 



expand all products, i.e., we obtain 
'H 6 (l-x 



M 



c + sx 



/ t—' \c + sx 

' 1=1 v 



x in Hft(l — x) and 
(11.6) 



+ e 



,'H (x) 

Finally, perform the Mellin transforms M( „ g '„^ , n 



c+sx 



i.e., 



with e,ei G 

express it in terms of S-sums and constants in terms of infinite S-sums. Note 
that these harmonic sums are less complicated (see below). 
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Replace M 



(x) 



c+sx 



,n) in (11.4) by the result of this process. 



For all harmonic sums that remain in (11.4) apply the above points. Since these harmonic 



sums are less complicated (see below) this process will terminate. 



Concerning termination the following remarks are in place. Since ^ 1 in (11.4), we know 
due to Lemma 7.5 that mi ^ 1 in ( 11.4). If not all mj ^ 1 in ( 11.4), we have to transform x — > 1— x 
in H mi! m2,...,mi(%) as described in Section 3.2 Due to Remark 3.5 we see that the single generalized 



polylogarithm at argument x with weight /, which will emerge, will not have trailing zeroes since 
mi ^ 1. Hence we can compute the asymptotic expansion of the Mellin transform containing this 
generalized polylogarithm using repeated integration by parts. We might be able to compute the 
asymptotic expansion of some (or all) of the other Mellin transforms by repeated integration by 
parts as well, but if we fail we know at least that due to Remark |7.5| the .S'-sums which will appear 
in (11.6) are less complicated (they are even of lower weight) than S ai)a2i ... )(Ife (6i, &2j • • • ,bk',n) of 



(11.4) and hence this algorithm will eventually terminate. 

Applying the described algorithm to S2,iQ, |;rt) produces the following asymptotic expan- 



sion: 



82,1(3, h n ) 



Si, 2 (i;,3;oo) + Si(|; 
126 12 1 



00 



125n 4 25n 3 
1 / 4074 
6~ 



5n 2 
366 



33- 



~S 2 (i;oo) 



~)Tl — 1 

n 5 
42 



S 2 (4;oo) 



1464 
6 " I ~625n 5 
2 n-i- - 



3125n 5 625n 4 125n 3 



6 

25^ 



?(•' 



1 

hn 



S 2 (i;oo) 



28 



+ S 1 (i;oo) S 2 (i;oo) +S 3 (i;oo) + i - — 4 



1 / 2037 183 
6" V3125n 5 ~ 625n 4 
183 21 3 



21 



625n 4 

1 
5n 2 



125n 3 
log(2), 



125n 3 
25n2 ~ 10 



25n 2 



1 

I / 2037 



Wn) ^ + 6" ^ :VI.2.-»// r - 
1464 126 12 



625n 5 



125n 4 25r7, 3 



with S m (|,oo) = Li m (I) , m,k£N + ,m>2. 
11.4 Further extensions 

The following lemma extends the presented algorithm to 5-sums S ai> ... )afe (&i, . . . , bk\ n) provided 
that if 1 61 1 > 1, we have that \bi\ < 1 and if b% = 1 then aj > 1 for 2 < i < k. In this case, we 



obtain the required form given in (11.4) except that < c < 1 and that the occurring S'-sums in 



T satisfy again the constrains from above. Then the following lemma yields a suitable integral 



representation of M 



2i0r) 



c+sx 



n 



and repeated integration by parts produces its asymptotic 



expansion. In addition the occurring S'-sums in T can be handled by recursion of the extended 
method. 



Lemma 11.4. Let H mijTO2i .., jTrafc (a;) be a generalized polylogarithm with G 
for 1 < i < k and < c < 1. Then 



\ (-1,1] U {0} 



M 



H 



mi,m2,...,mi 



(X) 



, n 



c — x 



1 I x (H mi) 



C ./ 



m2,.-.,m k 



(x) H mij?Tl2i . mfc (c)) 



dx 



c — X 
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+H mi)C;m2i ... imfc (c) + Hmi,m2,c,...,m k ( c ) + " " " + ^m 1 ,m 2 ,...,m k ,c(c) 

~l~H m2) ... : m k (c) Ho j?Tll _ c (l — c) + H m3i _ jmfe (c) Ho i?Tll _ Ci?Tl2 _ c (l c) 

H 1" H mfc (c) Ho, mi - c ,...,m fc _ 1 -c(l — c) + Ho, mi - c ,..., mfc - c (l — c) 

— 5'i i {i}( ? T')H miim2i ... imfe (c) . 



Proof. We have 



with 



where 



A [( H mi , m2 ,„„ mfc (x) \ = I" 1 ((f)"- 1)H 

C-X J J C- X 

1 (g) (Hmi,m2,...,m t (^) ~ H mi ,m2,...,m t (c)) ^ 
/O C — X 



/' 

JO 



I { X c ) n ^-m 1 ,m 2 ,...,m k {c) H mi ,m 2 ,...,m k (x) ^ 
C — X 

f 1 (f )"(Hini,Mi2,..,m t ( 3: ) ~ H miim2v .. ;m , fc (c)) 
7o c ~ x 



,(<0 I 
Jo 



I (l) n -l 

+ Hm 1 ,rn2,...,m k (c) J '- — - dx 



L 



^jm 1 ,m 2 ,...,m k (x) ^m 1 ,m 2 ,...,m k {c) ^ 

o c-a; 



^4 — -Si (-; n^j H mi ,m 2 ,...,m k (c) 

g f ^rn 1 ,m 2 ,...,rn k {%) ~ ~H-rn 1 ,rn 2 ,...,m k (c) , 

Jo c-x 

V v ' 

Bi:= 



i "mi,m 2 ,..,m li M H miiTn2[ ,.. |mfc (c) ^ 
C — X 



B 2 := 



^1 — ^rni,c,m 2 , ... y m k (c) *** H rniJn2) ... : m k ,cip 
»l-c 



g _ I ~Rm l ,m 2 ,...,m k {x + c) H mi ^ ,.,,,m fc (c) ^ 

— X 



J 

Jo 

H m2i ... jmfc (c) H mi _ c (x) + • • • + H mi _ Ci . _ imfc _ c (x) 



/ 

■/ o 



= H m2ji ,, : rn k (.C) Ho 5 mi— c(l — c) ^im^,...,rn k {c) H() jmi _ Cim2 _ c (l — c) 
— • - - — H mfe (c) Ho, mi -c,...,m/ i; _i-c(l ~~ c ) ~~ Ho imi - c ,...,m fc -c(l ~~ c ) j 

thus the lemma follows. 

A typical example for this extension is, e.g., 

, x , 1331 63 23 5 ; 1 x , 
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S2 (|;o°) i 



with Si(i,oo) = ln(m) - ln(m — 1), m € N\{1}. 

More generally, this lemma enables one to deal also with S'-sums which are outside of the 
5-sum case. To be more precise, our strategy works if the given S-sum S air .. )afc (&i, . . . , b^; n) 

can be written in the form (11.4) where M^ Hmi '^f' gx ' m '^ , n) can be expanded directly by 



Lemma |11.4 In addition we need that T consists of 5*-sums which are less complicated than 
Sai ...,a k (bi, . . . , bk', n) and which can be handled by iterative application of our expansion strat- 
egy (and using Lemma [ll.4| if needed). 



E.g., our expansion technique produces 

Si.u f-i,2,i;n) ~Si,i,i'-i «> 1- 

2 



'1,1,1 

+ 



1 

2 



l27n 3 



_ 9 _ 
21 1 4 ' 

1 



) + (" 



1 

' 27n- 



+ 



1 
9u- 



6n) ( ! 

2 2 1 

27ra 3 + 9rf 



+ 



+ 



H) 
15n 3 



+ 



\3n 2 



27n 3 
2 

~ 9W 3 



+ 



9n 2 



+ 



1 

3n 

1 

2 



2 
1 

27n 3 



log 3 (2) 



1 

9r? 



+ 



+ ( 27n 3 



2 1 

9n 2 3n 



log(2)log(3) + (-l) 



(-1) 
V9n 3 



3/i 



1 

3n 2 



l s 2 

log(3) 



Si(2;n) 



log 2 (2) 



log(3). 



12 Application from Quantum Chromodynamics 



In the following we concentrate on the discussion of the emergence of generalized harmonic sums 
in calculations of massless and massive Wilson coefficients in deep-inelastic scattering in QCD. 
Up to the 2-loop level neither for the anomalous dimensions nor for the Wilson coefficients 
generalized harmonic sums occur, cf. [S3 E21 123 ES]- They are also absent in case of the 3-loop 
anomalous dimensions |91| 1126] . In intermediate results of the calculation of the massless 3-loop 
Wilson coefficients they emerge, e.g. for a graph called LA27box(iV) in Ref. [125]. The generalized 
sums (2.1) have numerator-weights over the alphabet X{ £ {1,-1,2,-2} up to weight w = 6. 
Examples are: 

Si(2;n), 81,1,1,1,2(2, 1, 1, 1, l;n), S 2 ,i,3(2, 1, -1; n), S 3 , 2 ,i(-2, -1, l;n), etc. 

All of these sums contain besides the numerator weights it £ {—1, 1} only either the weight 2 or 
—2 once. For these cases one always may represent the corresponding sums such that 

-1 (±2x) n 
ax 

'0 



Sni,. 



(±2,ii,.. 



f 1 dx i±2x)n 1 Yc j H d ..(x), 
J ±2x-l ^ 3 A h 



(12.1) 



where Ha-\x) are usual harmonic polylogarithms. Since the corresponding expressions are 
weighted by a factor 2 _n , the individual expressions are transforms by the non-singular ker- 
nels ((±l) ri x n — l/2 n )/(±2x — 1). It even turns out that all generalized sums cancel in the final 
result, |125| . In this way, the massless Wilson coefficients up to 3-loop order do not depend on 
generalized harmonic sums. 

In Ref. [3] the 0(n f T^C A ,F) contribut ions to the 3-loop massive Wilson coefficient contribut- 
ing to the structure function F2(x,Q 2 ) at high virtualities Q 2 3> m 2 were calculated. Here m 
denotes the heavy quark mass. In this calculation most of the generated nested sums were calcu- 
lated directly using the package Sigma (1 10|, 11121 1 1111 1113| 1115] . In this way individual Feynman 
diagrams were split into many terms being summed individually. Here the contributing alphabet 
of numerator weights was Xi E {1,— l,j,— i, 2, —2}. Examples of generalized sums are : 
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Si (|;n) ,S 2 (-2;n),S2,i(-l,2;n),S 3j i (-2, -\\n) ,81,1,1,2 (-1, |,2,-l;n) , S 2 ,3 (-2, -\\n) , 
S2,2,i (-1, -5, 2; n) , 81,1,1,1,1 (2, 5,1,1, 1; n) , etc. 

The highest contributing weight was w = 5. In the final result all generalized harmonic sums 
canceled for the 0(njTj,C J 4 i i?)-terms, even on a diagram-by-diagram basis. 

This brought up the idea to combine all sums contributing to each diagram leading to 
C. Schneider's package SumProduction [119] . Through this a much smaller amount of nested 
sums, however, with very large summands, have to be calculated, [301 E]. We have applied this 
method for the gluonic 0(rifTpCA,F) contributions in [8]. 

In case of the massive 3-loop operator matrix elements generalized harmonic sums do not 
vanish always. For a series of ladder diagrams with six massive propagators they contribute for 
scalar integrals up to weight w = 4, [7]. This is likewise also the case for a series of Benz-type 
diagrams, [9]. Examples are: 

Si (2;n) ,Si, 2 (5, l;n) ,Si,i,i (|,1, l;n) ,Si,i, 2 (|,1, l;n) ,81,1,2 (2, 5, l;n),Si,i,i,i (2, §, 1, l;n) ,etc. 

Representing results for individual Feynman diagrams containing also 5-sums sometimes involves 
divergent S-sum expressions in individual terms as n — > 00. On the other hand, combinations 
of different terms contributing to the corresponding diagram are always convergent, as also the 
whole diagram, see [TtllOlfTT]. We therefore limit the consideration to convergent combinations of 
terms containing 5-sums here. Let us consider an example of Ref. [7\. Integral Ii(n), Eq. (3.18), 
contains the following terms growing like 2 n for n — > 00 : 

m , % r, n 3 — An 2 — 25n — 28) r „ /1 \ „ ,-1 \ 1 

Tl{n > ~ 2 (u + !)*(„ + 2) (n + 3)3 ^ «• ^-O + ^ «• " 2& 1 

T ' An) = ~\ (n + X)(l+3) I 5 '' 1 ' 1 ' 1 ( 2 <h 1 < 1 '- n ) +2S U-2 (2,l,l;n) -2C 3 S 1 (2;n)] . 

One may represent these sums in terms of Mellin transforms, cf. Eq. (3.36)-(3.40), [7]. There it 
turns out that each of the S-sums contributing to T\(n) can be expressed in the form 

n,i(3 + ^Fi,i(n), 

where n j G Q and F\^{n) are convergent integrals for n — > 00. Replacing the 5-sums with these 
alternative representations in T\(n) shows that the £3 part vanishes. In addition, rewriting the 
introduced integral representations back to 5-sums produces finite S*-sums. Similarly, each of the 
iS-sums in T2(n) may be expressed in the form 

r2,iC3<S , i(2;n) + F 2> i(n), 

where r2,i G Q and -F 2 ,i(n) are convergent integrals. Again, replacing the S*-sums with these 
representations in T2 cancels the S*i(2;n)^3 part; in addition, the introduced integrals can be 
transformed to an expression in terms of absolutely convergent S-sums. The remainder terms of 
l4(n) are even nested harmonic sums. Due to this only absolutely convergent S'-sum contribute 
to 14(11) via T\(n) and T\(n). A decomposition of this or a similar kind is expected in case of 
physical expressions in general. As the above example shows a suitable decomposition of the 
contributing S'-sums may be necessary. 

In Refs. [101 E] we a l so considered graphs with an internal 4-leg local operator insertion and 
five massive lines. Here generalized harmonic sums of up to weight w = 5 occur. Examples are : 

Si, 4 (|, -2; n) , Si, 2 ,2 (5, -2, -1; n) , Si, 2 ,2 (5, -2, 1; n) , Si, 2 ,2 (|, -2, -1; n) , etc. 

For this diagram also generalized harmonic sums nested with binomial- and inverse-binomial 
expressions contribute. 
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Generalized harmonic sums are a necessary asset both in higher loop and multi-leg Feynman 
diagram calculations. They are connected to the generalized harmonic polylogarithms via a 
Mellin transform. We have worked out their various relations and considered special new numbers 
associated to these quantities in the limit n — > oo of the generalized harmonic sums and special 
arguments of the generalized harmonic polylogarithms, generalizing the multiple zeta values. 
The relations discussed generate basis representations in the respective algebras. Note also that 
special classes of higher transcendental functions are contained in the generalized harmonic sums 
for n —7- oo. 

In cases that generalized harmonic sums occur in physical results an efficient way of imple- 
mentation consists in working in Mellin space, as being the case for the nested harmonic sums (see 



also the comments on page 31 ). Having determined the singularity structure of the problem, and 
knowing the shift relations n — > n + 1 of the corresponding quantities in analytic form, the knowl- 
edge of the asymptotic expansion at \n\ 3> 1, n E C is sufficient for the analytic continuation of 
anomalous dimensions and Wilson coefficients, similar to the case of the nested harmonic sums. 
The asymptotic representation is requested to be free of exponential divergences. Logarithmic 
terms ~ ln fc (n) with k related to the number of loops may contribute, however. The structure of 
the singularities may become more general, however, and has to be determined for the respective 
problem under investigation. A single numerical contour integral around the singularities yields 
the inverse Merlin-transform from n— to x-space. cf. [32J. The techniques outlined above also 
allow to perform the inverse Mellin-transform analytically for applications in x-space directly. 
For this purpose the corresponding basis-functions still need to be represented numerically at 
sufficient precision, unlike the case in n-space. 

13 Appendix: Available commands of the package Harmonic- 
Sums, m 

In the following we list some of the main commands of the code Harmonic Sums, a pack- 
age implemented in Mathematica, and illustrate their use by examples, with emphasis on the 
case of S'-sums. The HarmonicSums .m package and the conditions of use are found under 
www. rise . jku. at/research/combinat/sof tware/HarmonicSums/. 

S[ai, . . . ,a m ,n\ : defines the harmonic sum S aij ... jam (n) of depth m with upper summation 
limit n and a* E N, while CS [{{ai, b±, c{\ , . . . , {a m , b m , c m }} ,{si , . . . , s m } ,ri] defines cyclotomic 
harmonic sums S( ai fii,ci), ...(a m ,fem,cm)( ,Sl ' • ■ • > Sm > n 

) with Oi, Cj E N, hi E N U {0} and s; L E {-1, 1}; 

see p2 (2.1)]. 

S [ai , . . . ,a m , {xi , . . . ,x m y ,ri\ : defines the 5-sum S aij ... )am (xi, . . . , x m ;n) of depth m with up- 
per summation limit n ancj^jaj E N, Xi E M*. 

Z [ai , . . . , a m , {xi , . . . , x m y , n] : defines the Z-sum Z ai a m (•£]-) • • • jX m ]n) of depth in with 



upper summation limit n and «j E N, xi E M*; see (9.1). 

H[mi, . . . ,m n ,x'\ : defines the harmonic polylogarithm H mi) ... jmn (x) (mi E {—1, 0, 1}) or a gener- 
alized polylogarithm H mi ,...,m„( x ) ( m i ^ while H [{ai, &i}, {a n , b n },x] defines a cyclotomic 
harmonic polylogarithm [5] with di,bi E N. 

SToZ [S la% , . . . ,a m , {xi , . . . ,x m y,ri] ] : converts the S-sum S ai) ,„ >am (xi, . . . , x m ; n) to Z-sums. 
Example: 

SToZ[S[l,l,{l/2,-l/4},n] = Z 2 (~,nj + Z hl (l,-\,n 



Obviously, the Xi must be chosen from a subfield of R which is computable, in particular which can be treated 
within the computer algebra system Mathematica. 
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ZToS [Z [ai , . . . ,a m , {xi , . . . ,x m y ,ri] ] : converts the Z-sum Z 0lj ... j0m (x±, . . . , x m ; n) to S'-sums. 
Example: 

ZToS[Z[l,l,{l/2,-l/4},n] = -S 2 + S M Q,-^," 



ReduceToBasis [expr] : reduces occurring harmonic sums, S-sums or cyclotomic sums to basis 
sums which are (in the setting of their quasi-shume algebra) algebraic independent over each other. 
In general, tables (available at the HarmonicSums homepage) are exploited for this reduction. If 
sums are not covered by the given tables, one can set the option Dynamic — > Automatic: then 
for such sums (outside of the tables) the reduction is calculated online. Using the setting Dynamic 
— > True the reduction for all sums is calculated online. 

Example: 

ReduceToBasis[S[l, 2, 2, n], n, Dynamic — > True] = 

S 2 (n) (S 3 (n) - S 2>1 (n)) - S 2 , 3 (n) - 5 4 ,i(n) + S 2) 2,i(n) + *Si(n) (S 2 (n) 2 + S 4 (n)) + S 5 (n) 



ReduceToHBasis [expr] : reduces occurring harmonic poly logarithms, generalized polylogarithms 
or cyclotomic polylogarithms to basis polylogarithms which are algebraic independent over each 
other by using given tables. If the option Dynamic — > True is set, the reduction is calculated 
online; this will help if certain polylogarithms are not stored in the table. 

Example: 

ReduceToHBasis [H[l, 0, -1, l,x]] = 

(H (x)H 1 (x) - H 0A (x)) + H- ltl (x)H 0tl (x) - H^x) (H^(x)H 0:1 (x) - H - ltl (x) 
-Ho^-iix)) - H^x) (H (x)H_ 1A (x) - H^(x)H 0A (x) + tf ,i,-i (x)) - 2F ,_ 1 , 1)1 (x) 
-H 0>1 -i tl (x) 



Transf ormToSSums [expr] : returns an expression which transforms indefinite nested sums to 
harmonic sums, S'-sums or cyclotomic sums whenever possible. 

Example: 



expr = (n - l) 2 £ ^ + 4 C s(n - 2)\n - if ± U£ + n\n + if ± 



i=4 



i=l 



i=2 



^ 1 ^ (-3V 3, . 2 , , 3 A 1 ^ " (_i)< 

+£^E(iVw-2>" 2)(1+n) ^ 

1=1 j=l v ' i= 4 j=4 v «=1 



+E 



(-3)* 



-(l + 3j) 2 2 
rSi(*) 



Transf ormToSSums [expr 
1 



5(i A i),(3,i,2)(l, -3;n) + 5(3 ) i i 2),(i ) o,i)(-3, l;n) + 2(n - 2) 2 (n + l) d S_ 3 (oo) x 

(-1)" (-1)" , c , , , S_ 2 (n) , c , , 00 , , , (-1)"- 1 , 2(-l)' 



-Si _ 2 (n) - 3 



+ S_ 3 (n) + 



n 



+ S_ 2 (n) -2S- 1 (n) + 



n — 1 



+ 



n 



(-1) 



n— 1 



(n- l) 5 



+ 5_ 2 ,i(n) + (n-2) 2 



27 



:S 2 



1 \ (-l)"3" n - 3 (-l)™"^-"- 1 
3 ;n J + ^ + (n - 2) 2 
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n— 1q— n— 2 



+ 



(n-1) 2 
(-If 



+ n 2 (n + 1) E 



■1)' 



fiL 3 (n) + 4(n-2) 2 (n-l) 5 



5-i (n) 



-1) 



n-l 



n — 1 



n 



^(oo) 



Here S(ai,6i,ci),...,(a ra ,b m ,c m )( a; i> •••> x m, n) are generalized cyclotomic sums, cf. Ref. [S]. 

ReduceSums [expr] : returns an expression which transforms indefinite nested sums to harmonic 
sums, S'-sums or cyclotomic sums whenever possible and reduces occurring harmonic sums, S- 
sum or cyclotomic sums to basis sums which are algebraic independent over each other by using 
given tables. If the option Dynamic — > True is set, the reduction is calculated online; this will 
help if certain sums are not stored in the table. 

Example: 



ReduceSums [expr] 



1 



^ S 5(i,o,i),(3,l,2)(l) -3; n) + 25(3^2), (i, ,i)(- 3 > 1; ") - 3(n - 2) (n + 1) S 3 (oo) 



(-l) n (-l) n 9<? r s/ s/u l ,\ (-1)" , 2(-l)" , ("l) n 
— 5 2 26_i(nj + S_ 2 (n) I — Si(n) + — + 1 ) H h 



71° 



+ 



n — 1 n (n — l) 2 
2 



3-"- 3 (2(-l) n (7n 4 - 12n 3 + 10n 2 - Yin + 4) - 2 3 n n 2 (n 2 - 3n + 2) 2 S 2 (-£; 



+2n 2 (n+l) 3 ' ( ^ 



5-s (n~ 



71" 



(n — l) 2 n 2 

(n - 2) 2 (n - 1) ((n - l)w5_i(w) + (-1)") S 3 (oo) 

n 



+2S_ 2jl (n)j 

LinearHExpand [expr] : expands all the products of harmonic polylogarithms, generalized poly- 
logarithms and cyclotomic polylogarithms at same arguments. 

Example: 

LinearHExpand[H[0, l,x] * H[-l, 0, -1/2, x]] = 

ff-i,o,-!,o,iO) + 2-^-1,0,0 -l.i^) + 2jff -i,o,o,i ,-§0*0 + H -1,0,1,0,- 1 ^) + ^0,-1,0-1,1^) 
+-^0,-1,0,1,- 1 ( x ) + #0,-1,1,0,- + #0,1,-1,0,-^) 



LinearExpand [expr] : expands all the products of harmonic sums, S'-sums and cyclotomic sums 
at same arguments. 

Example: 



LinearExpand[S[l,2,n] * S[-l,2, 1/2, 1/3, n]] 



1 1 



2'3 



1 1 



5o,4 7T, , ; n) - 5-1,1,4 o' o ' n ~ -5-1,3,2 o' q ' 1 ! n ~ ^o^ 2 



2' '3 



1 1 



2'3' 



1 1 

2'3 : 



1; n 



-5 ,2,2 Q, 1, - 5i _i )4 ^1, i, - Sx tl)2 (l, ^, + 5.1,1,2,2 Q, 1, ^, l;n 

+5-i,i,2,2 Q, 1, 1, rt"j + 5-i,2,i,2 Q, 1, 1; n) + 5i -1,2,2 f 1, ^> ^> 1 : " 
+5i -1,2,2 f 1 ' 2' 1 ' 3 ;n ) + Sl > 2 - 1 ' 2 ( 1 ' 1 '2'3' n 
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InvMellin [expr, n, x] : calculates the inverse Mellin- Transform of expr. 
Example: 



ReduceConstants[lnvMellin[S[2, 2, 1,{1, 1, l/2},n],n,x],ToKnownConstants — > True] 

2-«5 2>2 (|,l;oo) 2-"- 1 (C2-ln(2) 2 )if , 2 (x) ln(2)2-»flo,2,o(z) , ln(2)tf ,i,o(x) 



x-2 x-2 x-2 

2- g Q,,,, 2 (x) 2-5 4 (l,co) ffig - ^ + ¥ 



+ 



x-2 



x-2 



x-2 



+ 



x — 1 
8(x- 1) 



x — 1 x — 2 



^ 5 2 " ^V/ijixi ln;2; - ^ 6 2 1 s ) 



x-2 



+S 4 -;oo +ln(2) 



x — 1 

/ln(2) 3 ln(2)C 2 7( 3 
1 6 2 8 



1 



2 - x 
^2^2 



ln(2)(-(C 2 -ln(2) 2 ) z -2(-5 1 , 3 (-,l;oo 



4 

Sl,4 



1 



l;oo - S; 



3,2 



-, l;oo 



+ ■^..2.2 ( ^,l,l;ooj + S 5 ( -;oo 



For the command ReduceConstants we refer to page 67 Note that some of the constants may 
be further identified by 



S k {z,oo) = Li fc (z) , \z\ < 1 . 

Mellin [expr, x, n] : calculates the Mellin- Transform of expr. 
Example: 



Mellin[H[2, 1, -3,x],x,n] /. S[l, {-1/3}, Infinity] -> -Log[4/3] = 

2" +1 5i, 2 Q,-§;oo) 5i, 2 (i-§;oo) _ 2" +1 5 2 ,i (l, -|;oo) S 2 ,i(l,4;oo) 

n+1 n+1 n+1 n+1 

2«+ 1 5i,i,i(i,2 J -i;oo) 5i,i,i(i,2,-|;oo) 5i >x (-3, -§; n) 5 3 (-|;oo) 



+ 



n + 1 



n+1 



(n + 1) 1 



n + 1 



2" +1 5i,i,i (i, -3, -|; n) 2*+* log (|) (|, -3;n) 2^ log (|) 5i tl (|, 1; n) 



n + 1 



n + 1 



n+1 



| 2" +1 53(-^oo) + (-l)"3" +1 5i (4;n) _ log (f) 5i(-3; n) + log (|) S^n) 



n + 1 



(n + l) ; 



(n + l) ; 



(n + l) ; 



1 | (-lrS^logd) | log (|) 



(n+1) 



(n+ 1) J 



(n + iy 



Dif f erentiateSSum [expr, n] : differentiates expr w.r.t. n. 
Example: 



DifferentiateSSum[S[4, 1, {l, -l/2},n],n] 



1 



1 



1, -x; n - Sa 2 1, --; n - 455,1 1, --; n + #-1,0 



f ^0,0,0,1,0(1) + #-1,0,0,0,-1,0 



1 



S4 (n) 



The constants may be further reduced to known multiple zeta values and/or cyclotomic numbers 
applying the command ReduceConstants [expr] (see page 67). 
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RemoveLeadingl [H [mi , . . . ,m n) x]] : extracts leading ones of H mii ... ;mn (x). 
Example: 

RemoveLeadingl[H[l, 1,-1,0, 1/2, x]] 

= ^iW 2 ff 1A iW - H^H^i^x) - H^H^i^x) - H 1 (x)H_ 1 ^i(x) 
+H_ 10 i 1 +H_ 1Q1 i ± (x) + H_ 1Q11 i(x) + H_ 110 i + H_ 1101 i(x) 

RemoveTrailingO [H [mi , . . . ,m n ,a;]] : extracts trailing zeroes of H mii ... jmn (x). 
Example: 



RemoveTrailingO[H[l/2, -1/3, -1, 0, 0, x]] 

= \h q {x) 2 Hi _i - tf (z)#o i _i _i(x) - H (x)Hi _i - flo(aOffi _i _i(x) 

2 2' 3' '2' 3' 2> 3> ' 2> > 3' 



+^ooi -i -i( x ) + H 0± -i o-iW + ff oio-i + #i - ioo-i( x ) + ^1 o-i 0-l( 

' ' 2 1 3 ' ' 2 5 3 ' ' ^ u ) 2 5 ' 3 ' 2 ' 3 ' ' 5 1 2 5 ' 3 ' ' 

+^i o-i -i( x ) 

2 > U 1 U 1 3 ) A 

RemoveLeadinglndex [H [mi , . . . , m n ,x]] : extracts the leading index mi of H mii ... imn (x) 
RemoveLeadinglndex [expr, 6] removes the leading index b of all polylogarithms in expr. 

Examples : 



■r 



RemoveLeadingIndex[H[l/2, 1/2, -1, 1, 1/3, x]] = 

-Hi(x) 2 H_ 1>lt i(x) - Hi(x)H_ h i x i(x) - Hi(x)H_ hl i i(x) - Hi(x)H_ lh ii(x 
+H_, 1 1 , 1 (x) + 1 1 1 1 (#) + -Hi 1 1 1 1 (#) + 1 1 1 1 (x) + .Hi a! 1(1) 

>2>2> '3 '2' '3'2 '2' >2>3 ' >3>2>2 > >2'3'2 

+H_, , 1 1 

* ' x » 2 » 2 » 3 



RemoveLeadingIndex[H[b, 1/2, —1, 1, — 1/3, x], b] = 

H i -1 1 - Hi _ ±1 _i b (x) - iii _! 16 _i(x) - #1 _ lbl _i(x) -Hi b _ 11 _i(x) 

2' ' ' 3 2' ' ' 3' 2' ) 11 3 2' 5 ' ' 3 2' ' ' ' 3 

RemoveTrailinglndex [H [mi , . . . ,m n ,x]] : extracts the trailing index m n of H mii ... )mn (x). 
RemoveTrailinglndex [expr, 6] removes the trailing index b of all polylogarithms in expr. 

Examples : 

RemoveTrailingIndex[H[-l, 1, 1/3, 1/2, 1/2, x]] = 

\H h (xfH_ lM (x) - H h (x)H_ 1>hl; >(x) - II Lnll ,,.(,•) - H^H^^x) 
+H_, i i , i (x) + H_, i i i i (x) + H_, , i i i (x) + Hi i i i (x) + Hi _, , i i (x) 

1 '2'2' '3 '2' '2'3 ' '2'2'3 2' '2' '3 2' ' >2'3 

+Hi i , i(x) 

2'2' ' '3 

RemoveTrailingIndex[H[-l, 1, 1/3, l/2,b,x],b] = 

H_ ± i i i (x)H b (x) - H_! 1 i b i(x)- H_ x lh ii[x)- H_ x b 1 i i (x) - H b _ x ± i i (x) 

' '3'2 ' '3' 1 2 1 ' '3'2 ' ' '3'2 ' ' '3'2 

SRemoveLeadingl [S [a\ , . . . ,a m , {x± , . . . ,x m } ,n] ] : extracts leading ones of 

Sai,...,a m (xi, . . . , x m ; n). 

Example: 

SRemoveLeadingl [S[l, 1, 1, -2,4, {1, 1/3, 1/2, 1/2, 1/2}, n]] = 
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5i(n)Si,i _2,4 Q, ^ n ) + -2,5 Q. 



1111 \ _ /l 1 1 1 N 



A i i i \ /i i i i \ /ill i 

+51,2,-2,4 (j, 3. 2' 2 5 "J + 52 ' 1 '- 2 ' 4 2' 2' 2 ;n J " 5l ' 1 '- 2 . 1 ' 4 ^3' 2' 2' 2 5 

-51,1,-2,4,1 (^^4 ^ ^ " 5l ' 1 ' 1 '- 2 ' 4 ^ " 

-51,1,1,-2,4 (3' 1 '2'2'2 ;n 

SRemoveTrailingl [S [ai , . . . ,a m , {xi , . . . ,x m } ,n\ ] : extracts trailing ones of 

^ai,...,a m { x li • • ■ > -^mi 

Example: 

SRemoveTrailingl[S[-2,4, 1, 1, 1, {1/3, 1/2, 1/2,-1/2, l},n]] = 
„ , , /111 1 \ /111 1 \ /111 1 

5i(n)S_ 2 ,4,i,i ^3, ^> 2' " 2 5 n ) + 5 ~ 2 ' 4 + 2 2' 2' " 2 ; "J + ^A 2 ' 1 \3' 2 ' 2' " 2~ ; " 

/111 1 \ /111 1 \ /111 1 

+5-2,5,1,1 ^3, ^ 2' " 2 ; n ) + S - 1 ' 4 ' 1 ' 1 U' 2' 2' " 2 ; "J " 5 " 2 ' 1 . 4 ' 1 ' 1 U' ' 2' 2' " 2 ; " 

-5_2,4,i,i,i Q, 1, -\, n) - 5-2,4,i,i,i Q, ^ 1, ^, - ^ n 

f I I I I 

-51,-2,4,1,1 ^'3>2'2'~2 ;n 

SRemoveTrailinglndex [S [ai , . . . ,a m ,{xi, . . . ,x m },n]] : extracts the trailing index (a m ,x m ) 
of S ai ,..., am (xi, . . . ,x m ;n). 
Examples : 

SRemoveTrailingIndex[S[ai, a m , {xi, x m }, n]] = 

5 aa (x 3 ; n) 5 aij02 (xi, x 2 ; n) + 5 01)a2+03 (xi, x 2 x 3 , n) + S ai + a3 ,a 2 (xix 3 , x 2 ; n) 

— 5 ai ,a3,a2 O^li -^3) 2-2 j ^) 5 a3 , ai , a2 (x 3 , Xi, X2! 7l) 

Transf ormH [H [mi , . . . ,m n ,y\ ,x] : performs several transforms on the argument of a general- 
ized polyloarithm. 

Examples : 

Transf ormH[H[-l, 2, 0, 1/x], x] = 
# i (1) + tf , 2 (l)) (-#_i(x) + # (x) + ln(2)) - #-1,0,0^) - #-i,i, o 0*0 + #o,o,o(s) 

+i? 0ii0 (x) + i?_i,i, (l) + #-1,2,0(1) - #o,|,o(l) + X 

Transf ormH[H[- 1,2,0, 1 - x],x] = -H 2fi (l)H 2 (x) - # 2 ,-i,i(x) + #-1,2,0 (1) 
Transf ormH[H[- 1,2,0, (l - x)/(l +x)],x] = 

#2,o(l) (-#-i(x)) + #-i,-i,-i(x) + #-i,_i,i(x) - #_ 1; _i ; _ 1 (x) - #- 1 ,_i,i(x) + #-1,2,0(1) 

HarmonicSumsSeries [expr, {x, c, ord}] : tries to compute the series expansion of expr with re- 
spect to x about c 6 1 U co up to order ord £ N. expr may contain harmonic sums, cyclotomic 
sums, and their generalizations as well as their associated poly logarithms. 

Examples : 

HarmonicSumsSeries[H[-2/5, 3/24, 0, 1/4, x], {x, 0, 10}] = 
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40x 3 4 2026a; 5 20081a; 6 12632239a; 7 1205544979a; 8 119322150247a; 9 

h 35x H 1 1 1 1 

3 9 18 1890 30240 476280 

188972192341a; 10 
+ - 



117600 

HarmonicSumsSeries[S[l, — 2, 1, n], {n, 0, 4}] = 

n fC^S _ 5|5 V n 2 (jig + CJ \ + n 3 (£ (oo) + 5„ 4 ,_ 3 (00) - 5_ 4 ,i,_ 2 (oo) 



168 4 

ClCs 29C 2 Cs >\ , 4 f q f x o r ^ 311n(2)C 2 C5 

- C7 I + n I S-6 -2(00) - 6_ 5 -3(00) + 6_ 5) i ,-2(00) - 



32 32 / \ ' K ' ' v 7 ' ' v 7 16 



x? , v , , 23C 2 4 3C2CI , 155C3C5 
+6 8 (oo) + s 6 C2 + ^ 7" + 



70 4 128 
HarmonicSumsSeries[S[l, — 2, 1, n], {n, Infinity, 10}] 



, 5 5 15 1 5 

+C3 [-7T- + T^-o - W7-I + T^TT"* " 777"* + 



2n 12ra 2 24n 4 252n 6 48n 8 132ra 10 

_ 3Cf , n n/ r _2L 1^11 59 118529 130213 \ „ , , n 

5 +( j V + « 4 3n5 12n6 + 60n 7 + 8n 8 630n9 2520n 10 ) 5Csin(nj 1 

HarmonicSumsSeries[S[l, -2, 1, {1/2, 1/3, -1}, n], {n, Infinity, 3}] = 

ls v (-1,2; 00) - ?S M (-1,6; 00) + | ta(2)S, (-i;oc) - ?ln(2)S, (-;U) 

+ 3 l„(2) Sl C 1 ; 00) - 92 -*H;°°) + 32 "' S 'H ; °°' - 32 ""' 5 '<^ ;C °) 
2 \6 / n 6 n z n 

3/1 \ / 1 \ 1011 / 1 \ 3 / 1 \ 71n(2)2-«3 2 -™ 
"tt! 00 ^1 _ 77;°° ^TT^i -77;°° _ 7762 — 7S00 - 



2 V 3 ' 7 V 6 ' 7 800 V 6 ' 7 2 V 3 ' 7 !25n 3 
ln(2)2 -„ 32 -„ _ 7 ln(2)2 _„_ l3l _„ _ 1 ln(2)2 -„- l3 -„„ _ ln(2)2-»-3-» 21 ln(2) 



25n 2 25 w 5 w ' 5n 32 

+ 79(-|)"2- n - 5 + 63 2-"- 4 21 2-"" 4 (_i)n 2 -n-4 3 i-n ^ j / l\ n 2 _ n _ 3+ (-|)"2 



343n 3 ra 3 ra 2 49n 2 7 \ 3 / 7n 

21 2-"-« _% _9_ /j, 
n 4 + 280 V 

with ln(6) = ln(n) + je- 

ComputeSSumBasis [w, {x±, X2, •••}, n] : computes a basis representation of 6-sums with indices 
Xj at weight w using algebraic, differential and half integer relations. 

Example: 

ComputeSSumBasis [2, {l, 1/2, 1/3}, n] = 

j j6_2(n),6 2 (n),6_i,i(n),6 2 Q; n ) > ^ Q ;n ) ,Sl ^ ( 1 '^ ;n ) ,Sl ' 1 ( 1 ' 2~' 
jfi-i _i(n) -> -5-1,1 (n) + 6_ 2 (n) + 5_i(n)6i(n), 
^i(n)^^6i(n) 2 + ^M, 
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( hn) ->• -5*1,1 U, + S^Si Q;™J + 52 Q;n j • 



Sj.-j ( -»• "5i,i +^i(n)Si Q;n) + 5 2 Q;n 

5 2 (n) 



S-i,-i(n)-^2S-i(n)^+ 

ComputeHLogBasis lw, x, Alphabet — > {ai, a^}] : computes a basis representation of general- 
ized polylogarithms with indices ctj at weight w using algebraic relations. 
Example: 

ComputeHLogBasis [2, x, Alphabet ->■ {1/2,-1/3,1}] = 

{{" .. : . j- 

j#i,i(x) *#i(x) 2 , 

if, i (x) -»• Hi (x)H 1 (x) - Hi , (x), 

1 '2 2 2 ' 

#i,_i(z) -> #_i(x)i?i(x) - ^_i 5l (x), 
#i,_i(x) -> iJ_i(x)i?!(x) - #_i,i(x), 

ifi i(x)^Jtfi(x) 2 , 

2 ' 2 Z 2 

fli _i(x) ^ H_i(x)Hi{x) - H_ii{x), 
H h _[{x) hJ(x)Hi(x) - H_li{x), 

Hi i(x)->-H i(x) 2 , 

ff-i,-i(x)-^^_i(x) 2 || 
ComputeSSumlnf Basis {xi, .., x&}] : 

computes relations of S'-sums at infinity and weight «; with indices in {xi, ...x^}. 
Example: 



5 , _ 2 (oo),S'2 ( -^;oo 



ComputeSSumInfBasis[2,{l/2, -1, -1/3, 1}] = 
1 

jsi,i (l,^;oo) -> 5i(oo)5i Q;oo) - ^-i(oc) 2 , 



3' J -^-'"i y 3 '~y 2 1 V 3' 

51,1(00) ->■ ^S'i(oo) 2 - S_ 2 (oo), 
Si,_i(oo) -> S_i(oo)5i(oo) - ^_i(oo) 2 , 
S 1 !,! f-,l;ooj ->• -S_ 2 (oo), 
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Si i I — , 1; oo i — > -S\ ( — ;oo| 2 + So ( — ; oo 
' V 3 J 2 V 3 / V 3 

5_i,i(oo) ->■ -SLi(oo) 2 + S_ 2 (oo), 
5_i ,_i(oo) -> -5_i(oo) 2 - 5_ 2 (oo), 
Si Q ;o °) -^ 5 '-i( 00 ) 2 ~ 5-2(00), 
5 2 (oo) -> -25_ 2 (oo)|| 

ReduceConstants [expr] : reduces (as much as possible) occurring harmonic sums, 5-sums or 
cyclotomic sums at infinity and harmonic polylogarithms, generalized polylogarithms and cyclo- 
tomic polylogarithms at constants to a set of basis constants using given tables. Via the option 
ToKnownConstants — > True a special set of basis constants is used, i.e., log(2), £ 2 , C3> • • • For the 
multiple zeta values we use the basis given before in |123| for weights up to w = 6. The relations 
were calculated newly and agree with those in |123| . Extended tables (for higher weights) both on 
the harmonic sums at N — > 00 and the the harmonic polylogarithms at x = 1 were given in asso- 
ciation with Refs. |123|, 11051 12H] and are available using FORM [124] codes. It is straightforwardly 
possible to generate more extensive tables for special numbers as harmonic sums, cyclotomic 
harmonic sums and classes of generalized harmonic sums at infinity and likewise special values 
of the corresponding harmonic polylogarithms at a given argument over respective bases using 
HarmonicSums. 

Example: 

ReduceConstants[S[4, 1, —1, Infinity], ToKnownConstants — > True] = 

Li (1} C + ln(2) ' C2 - 51n ( 2 ) 2 ^ 2 , 3m(2)C 2 C3 931n(2)C 5 _ 5se _ 149g 49Cf 
14 V 2 / 24 8 2 32 2 168 64 

Also tables of sums and polylogarithms at general argument are available. The harmonic sums and 
polylogarithms are tabulated up to w = 6, for the cyclotomic harmonic sums over the alphabet 
{±iy/i k , (±l) J /(2i + l) k } and the alphabet for the special generalized sums discussed in the 
present paper both to weight w = 4. Depending on the weight and length of the alphabet the 
corresponding calculations may become more demanding. The corresponding tables of relations 
grow correspondingly. If the option Dynamic — > True is set, the reduction is calculated online. 
This will help if certain sums are not stored in the tables yet. 

HToSinf [expr] : transforms generalized polylogarithms and cyclotomic polylogarithms at 1 to 
5-sums and cyclotomic sums at infinity. 

Example: 

HToSinf[H[l,2,-4,3, 1]] = 

1\„/12\„/14\„ / 1 2 



-Sifl [I, -; ooj - 5 2 , 2 -5 00 j - 5 3 ,i --; ooj + 8 W [I, -, -;oo 

+5i A1 (l -^-|;°°) + ^,1,1 Q (l. 2'"2'-3 ;0 °) +Sa 

Sinf ToH [expr] : transforms 5-sums and cyclotomic sums at infinity to generalized polyloga- 
rithms and cyclotomic polylogarithms at 1. 

Example: 

SinfToH[S[l, 2, 3, {1/2, -1/3, 1/4}, Infinity] = 
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#0,0,-6,0,0-24(1) — #0,0,0,0,0,-24(1) + #2,0,-6,0,0,-24(1) — #2,0,0,0,0,-24(1) 



HToS [expr~\ : computes the power series expansion of generalized polylogarithms and cyclotomic 
polylogarithms without trailing zeroes. 

Example: 

HToS[H[l,2, -4, 3,x] = 

°° 2-VS 2 (|. ri ) ~ 2 -^5i,i(-i,-|;Ti) A oc T ^S 3 (l;ri) 



n=l T ! n=i Tl n=i 



y- a; n 5i,2 (jj;n) A ^^2,1 (-j-|;n) _ A g^jW (i^iizilljj 
2^ r 2^ r 2^ r 

Tl = l 71 = 1 71=1 



" 4~^(-x)^5i (-fin) ~ 

2^ T 3 + 2^ r 4 

n=l 1 n=i 1 



SToH [expr] : takes (part of) a power series expansion of a (generalized/cyclotomic) poly logarithm 
and finds the (generalized/cyclotomic) polylogratithm(s) it originates form. It is the inverse of 
HToS. 

Example: 



SToH 



.71=1 



"Tl, 



3' 



-#[0, 0, -4, 3, x] + #[0, 0, 0, 3, x] 
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